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ABSTRACT 

Global stationary configurations of both aligned and logarithmic spiral magnetohydro- 
dynamic (MHD) perturbations are constructed analytically within an axisymmetric 
background of razor-thin scale-free gas disc, which is embedded in an axisymmetric 
gravitational potential of a dark matter halo and involves an isopedic magnetic field 
almost vertically threaded through the disc plane. The scale-free index /3 of the disc 
rotation speed vg oc i?"'' falls in the range of (—1/2 , 1/2) where R is the cylindri- 
cal radius. With the holding-back of a deep background dark matter halo potential, 
the isopedic magnetic field may be strong enough to allow for the magnetic tension 
force overtaking the disc self-gravity, which can significantly influence global station- 
ary MHD perturbation configurations and stability properties of the scale-free disc 
system. Only for stationary logarithmic spiral MHD perturbations with a perturba- 
tion scale- free index /3i = 1/4 or for aligned stationary MHD perturbations, can the 
MHD disc maintain a constant radial flux of angular momentum. The variable radial 
flux of angular momentum in the radial direction corresponds to a complex dispersion 
relation. The marginal instabilities for axisymmetric MHD disturbances are also ex- 
amined for a special case as an example. When the magnetic tension force overtakes 
the disc self-gravity, the scale-free disc can be completely stable against axisymmetric 
MHD disturbances of all wavelengths. We predict the possible existence of an isope- 
dically magnetized gas disc system in rotation primarily confined by a massive dark 
matter halo potential. 

Key words: galaxies: kinematics and dynamics — galaxies: spiral — galaxies: struc- 
ture — ISM: general — MHD — waves. 



1 INTRODUCTION 

As an important intermediate stage of many astrophysical processes, the dynamics of magnetized disc systems are frequently 
invoked to model different phenomena on various spatial and temporal scales, such as spiral or barred-spiral galaxies, magne- 
tized accretion discs surrounding supermassive black holes (SMBHs) anchored at galactic centres, accreting environs involving 
binary stars and protoplanetary systems etc. It is therefore of considerable interest to study magnetohydrodynamic (MHD) 
processes in discs for both theoretical and applied purposes. The pioneer development of the classic density wave theory for 
large-scale perturbation structures in differentially rotating discs (Lin & Shu 1964, 1966) has opened up a physical scenario to 
understand the basic dynamics of spiral galaxies (Lin 1987; Binney & Tremaine 1987; Berlin & Lin 1996). Linear perturbation 
or nonlinear theories are widely used as powerful tools of analysis to explore the structure and stability of disc systems. In many 
processes, initially small perturbations may lead to more violent consequences (e.g., gravitational collapse of a central core) 
which are crucial for the dynamical evolution of astrophysical systems. In various contexts, MHD perturbation configurations 
are valuable indicators of transitions to instabilities and of further dynamical structure and evolution. 

For simplicity and clarity, the scale-free disc system as an idealized model is a frequent target of theoretical investigation. 
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The concept of a 'scale-free' or a 'self-similar' disc means all relevant physical quantities in the disc system vary as powers 
of cylindrical radius R such that the disc system carries no characteristic length scale (e.g., the linear speed of disc rotation 
vg oc and the equilibrium disc surface mass density Eo oc R~^^~^ etc. with (3 being a constant exponent). Given possible 
situations, one may introduce 'boundary conditions' to describe a modified scale-free disc by cutting a central hole in a disc 
(e.g., Zang 1976; Evans & Read 1998a, b) or by limiting the radial extent of a disc. Observationally, many disc galaxies show 
more or less flat rotation curves, indicating the presence of massive dark matter halos according to the Newtonian gravity 
theory. Theoretically, discs with flat rotation curves (i.e., the azimuthal velocity vg = const or /3 = 0) belong to a specific 
subclass of scale- free discs, referred to as singular isothermal discs (SIDs). Since the seminal analysis of Mestel (1963), the SID 
model has attracted considerable attention in various astrophysical contexts (e.g., Zang 1976; Toomre 1977; Lemos, Kalnajs 
& Lynden-BeU 1991; Shu et al. 2000; Lou 2002; Lou & Fan 2002; Lou & Shen 2003; Shen & Lou 2004a,b; Lou & Zou 2004, 
2006; Shen, Liu & Lou; Lou & Wu 2005; Lou & Bai 2006). In most normal spiral galaxies, disc rotation curves tend to be 
flat over extended radii, giving a strong evidence for the existence of massive dark matter halos (e.g., Rubin et al. 1982; Kent 
1986, 1987, 1988). Given the inferred dominant masses, we recognize that a dark matter halo plays extremely important roles 
in the dynamical evolution of a disc galaxy (Miller et al. 1970; Ostriker & Peebles 1973; Hohl 1971; Miller 1978; Binney & 
Tremaine 1987). 

Besides the normal mode perturbation analysis in scale-free discs (e.g., Zang 1976; Binney & Tremaine 1987; Fan & Lou 
1997; Evans & Read 1998a,b; Goodman & Evans 1999; Shu et al. 2000), the zero-frequency neutral perturbation modes or 
stationary perturbation configurations are emphasized as the marginal instability in scale-free discs (e.g., Lemos et al. 1991; 
Syer & Tremaine 1996; Shu et al. 2000; Lou & Shen 2003; Shen & Lou 2003; Shen & Lou 2004a; Shen et al. 2005; Lou & Zou 
2004, 2006; Lou & Wu, 2005; Lou & Bai 2006). Axisymmetric instabilities are thought to set in through transitions across 
such neutral modes (Safronov 1960; Toomre 1964; Lynden-BeU & Ostriker 1967; Lemos et al. 1991; Shu et al. 2000). Parallel 
to zero-frequency modes for a^xisymmetric disturbances, Shu et al. (2000) tentatively proposed that logarithmic spiral modes 
of stationary perturbation configurations also signal onsets of non-axisymmetric instabilities; the resulting criterion appears 
to be compatible with the criterion of Goodman & Evans (1999) for instabilities in their normal mode approach. Recently, 
Shen & Lou (2004a) extended the analysis of Shu et al. (2000) in a gravitationally coupled composite scale-free disc system 
of gaseous and stellar discs using a two-fluid formalism. In solving nonlinear MHD equations, it is a challenge to achieve 
analytical results unless some simplified assumptions are introduced, such as the stationarity (uj = 0) condition etc. Most 
results are analytic in this paper and it is highly desirable that numerical simulations should be developed as touchstones for 
the analytical results and further venture into new regimes where analytical methods cannot reach (Syer & Tremaine 1996). 

Magnetic field receives considerable attention in astrophysics, and contributes to structures, dynamics and diagnostics 
on various scales of different astrophysical systems (Sofue et al. 1986; Beck et al. 1996; Balbus & Hawley 1998; Kaburaki 
2000, 2001; Balbus 2003; Vallee 2004; Hu & Lou 2004; Yu & Lou 2006; Yu, Lou, Bian & Wu 2006). Prompted by numerical 
simulations of magnetized cloud core formation through the ambipolar diffusion (e.g., Nakano 1979; Lizano & Shu 1989), Shu 
& Li (1997) introduced the so-called 'isopedic' condition requiring that the a— component flux {Bz in cylindrical coordinates) 
of an open magnetic field, which threads through the disc almost vertically, is proportional to the disc surface mass density 
(E) or Sz/E = const. With this simple and justifiable assumption (Lou & Wu 2005), effects of such an isopedic magnetic field 
can be simply subsumed into relevant terms in MHD equations (Shu & Li 1997); this simplification is powerful in analytical 
analysis in comparison with the usual awesome MHD equations describing the roles of magnetic field. Lou & Wu (2005) 
discussed the global MHD perturbation configurations in a composite SID system of coupled gas and stellar components with 
an isopedic (vertical) magnetic field. Lou & Zou (2004, 2006) investigated the similar kind of composite MHD SID problems 
except that their magnetic field is coplanax rather than isopedic. Lou & Bai (2006) analyzed MHD perturbation structures in 
a composite system of two coupled scale-free discs with a coplanar magnetic field. 

The purpose of this paper is to couple an isopedic magnetic field with a single scale-free gas disc embedded in an axisym- 
metric dark matter halo potential and to construct physically allowed global stationary MHD perturbation configurations in 
such a disc as well as to investigate the instability of such a disc system. Parallel to the hydrodynamic treatment of Syer & 
Tremaine (1996), our analysis leads to further MHD extensions and generalizations. To be specific, six components of forces 
are involved in the disc equilibrium of our analysis, namely, the radial centrifugal force, the thermal gas pressure force in 
the disc, the self-gravity of the disc, the holding-back gravity force of an axisymmetric dark matter halo and the tension and 
pressure forces resulting from an isopedic magnetic field within the disc plane. While highly idealized, the unique combination 
of these essential elements of a disc galaxy does provide a new perspective and has already shown some novel characteristics. 
For example, the introduction of an isopedic magnetic field geometry in a scale- free disc docs bring several new and interesting 
features in reference to prior related disc model analyses (Syer & Tremaine 1996; Shu et al. 2002; Shen & Lou 2004b; Lou & 
Zou 2004, 2006; Shen, Liu & Lou 2005; Lou & Wu 2005; Lou & Bai 2006). 

This paper is organized as follows. In Section 2, we describe the basic MHD model of a scale- free razor-thin disc. In Section 
3, isopedic magnetic field are discussed under the scale-free condition. In Section 4, we construct the stationary background 
axisymmetric equilibrium in a rotational MHD balance. In Section 5, a generalized MHD dispersion relationship is derived 
for MHD perturbations. We investigate the global stationary aligned MHD cases in Section 6 and stationary unaligned MHD 
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(logarithmic spiral) cases in Section 7. We summarize the results and provide discussions in Section 8. Specific details of 
derivations are contained in Appendices. 



2 A CLASS OF SCALE-FREE RAZOR-THIN DISCS 

Formally, a razor-thin disc has a three-dimensional mass density mathematically described by 

p{R, e, z, t) = E(i?, e, t)5{z) , (1) 

where cylindrical coordinates {R, 9, z) are adopted, 0, t) is the two-dimensional surface mass density and 5(z) is the 

Dirac 5- function with the vertical coordinate z as the argument. Except for the three-dimensional gravitational potential, any 
physical quantity 5 is a function of R, 9 and t as q{R, 9, t) for a razor-thin disc. 

In a razor-thin disc, the vertically integrated barotropic equation of state has a 'two-dimensional' form of 

n = feE^ , (2) 

where n(-R, 9, t) represents the two-dimensional gas pressure, constant coefficient k ^ and the barotropic index F > 0. We 
here take the general 'three-dimensional' polytropic equation of state p = kp^ for a comparison 



n(_R, 9, t)= / p{R, 9, z, t)dz ^ T,{R, 9, t)= / p{R, 9, z, t)dz , 

Jaz J Az 

where Yl is determined by vertically integrating the thermal gas pressure p. Given certain assumptions such as 11 oc 'Epc/ pc 
with pc and pc being the central pressure and density, a relationship between the two indices F and 7 can be established as 
F = 3 — 2/7 (see equation 3.5 of Lemos et al. 1991). 

The sound speed Vs in a barotropic disc can be defined as 

o dn 



and the enthalpy H is expressed by 



= fcFE^"^ = Fn/E (3) 

QZj 

^^-/- = ^. (4) 

y E (F-i) ^ ' 

The gravitational potential ^(R, 9, t) in the disc plane and the disc surface mass density is related by the Poisson integral 

$(i?, O^t)- j J [,2 + ^2_2i?,cos(e-,^)]i/2 ' 
where G is the gravitational constant. 

Among various model problems in disc dynamics, the scale-free disc model is often picked up by theorists for its relative 
simplicity and is explored as a powerful vehicle for a global analytical analysis. Scale-free discs carry no characteristic spatial 
scales. For a mathematical approach to scale-free discs, we adopt the definition given by Lynden-Bell & Lemos (1999). A fiat 
disc is scale-free if any physical quantity q{R, 9) of a disc configuration has the following relationship 

qiCR, 9 + h{(:)) = A^{0q{R, 9) , (6) 
where is a parameter, and h{(^) and Ai{(,) are two arbitrary functions. By requirement q has a general solution form of 

q = R°'F(9 + filnR) , (7) 

where a is complex, ^ is real, and F is an arbitrary function. The proof can be found in Lynden-Bell & Lemos (1999). 

Fluid equations provide a valid dynamical description for a gas disc component, while distribution function equations 
are more appropriate for a coUisionless stellar disc component (e.g., Binney & Tremaine 1987). After an introduction of the 
velocity dispersion of the stellar disc component regarded as an effective 'sound speed', ffuid equations can be approximately 
applied to stellar disc component as well (Lin & Shu 1964, 1966; Lin 1987; Binney & Tremaine 1987; Berlin & Lin 1996; Lou 
& Shen 2003; Shen & Lou 2004a; Lou & Zou 2004, 2006; Lou & Wu 2005; Lou & Bai 2006). 

In treating a razor-thin scale-free disc system, we adopt two-dimensional fluid equations. The three basic fluid equations 
of describing razor-thin disc dynamics can be written as 

fIT 1 r) ^ r) 

a^ + Sa^(^^-) + ^a^(^-)-°' 

dvR dvR ve dvR vj _ 1 dU c>($ + $) _ d{H + $ -f $) 

dt ^'"^ dR ^ R 09 R^ ■SdR dR ^ OR ' ^ ' 

dv0 dvg Vg dvg VgVR _ 1 dU 9(<E> + $) _ 9(_ff + <!> + $) 

"aF "^"^aR "R "90" R ^~T^~d9 Rd9 Rd9 ' ^ ' 

where vr is the radial bulk flow velocity and vg is the azimuthal bulk flow velocity and <1> represents the axisymmetric 



4 Y. Wu & Y.-Q. Lou 



background halo potential presumed to be unperturbed, i.e., P parameter remains unchanged in equation IIH . The last 
equalities in equations @ and llOII corresponds to the barotropic condition. 

In order to meet the scale-free requirement (Syer & Tremaine 1996; Goodman & Evans 1999; Shen & Lou 2004a; Shen, 
Liu & Lou 2005), we presume relevant physical quantities in the following forms of 

E = R-'^''^-^S{lp) , VR = R-'^a{(p) , V0 = R~'^b{(p) , 

$ = -i?-2'^P((^) , $ = -i?-2''P, H = R-^'^Q{ip) , (11) 

where we define ip{R, 9) = 9 + fj,\n R and /3 is a scale-free index. Both /i and l3 are constant real coefScients here. 
To satisfy the scale-free condition in equation we should require 

(-2/3- l)(r- 1) = -2/3 ^ r = (l + 4/3)/(l + 2/3) . (12) 

The constraint F > for a warm disc (fc > 0) implies either /? > —1/4 or 13 < —1/2 and thus ensures v'^ > 0. For a cold disc 
(fc = 0), such a constraint is unnecessary. 

Other constraints on (3 come from the disc mass distribution. First, E is assumed to decrease with increasing R and this 
implies /? > — 1/2 by the first relation in equation (IIH . Secondly, the central mass should be finite. This simply means 



2tt rr r^TT /■27r 



lim / / E(P,e»)7?dJ?de» = lim / / R'^'' S{9 + filn R)dRd9 ^ S(e')d6i lim 
>— 'OJo Jo ''^oJo Jo Jo >--.o 1-2/3 



< +00 





We then require 1 — 2/3 > 0, indicating /3 < 1/2. In short, we therefore have —1/2 < /3 < 1/2 based on the requirement of the 
disc mass distribution. 

Moreover, the allowed range of (3 is constrained by the convergence of the Poisson integral relating E and "1>. According 
to the analytical results of Qian (1992), for E(_R, 9) = /^-a/^-i+i"- exp(im6') = 71"^''"^ exp[i(me -|- i/\nR)] with u being a real 
constant parameter, we have correspondingly 

$(7?, 9) = -GYm {13 - ii^/2) _R-2''+i'^ exp(im6i) , (13) 
where Ym{x) is defined in terms of the standard F function as 

^F(m/2 - a. + l/2)r(m/2 + a:) 
' - F(m/2 - a; l)r(m/2 + x + 1/2) ' ^ > 

The above expression remains valid only for — m/2 < /3 < (m + l)/2; outside this range of /3, the force from gas materials at 
either r — > oo or r — > diverges (Qian 1992). When m ^ 1, potential expression l|18|l remains valid for /3 G (—1/2, 1/2). For 
the m = case, although the potential $ in equation H13|l converges only for /3 £ (0, 1/2), we can extend the valid range to 
/3 £ (—1/2, 1/2) if we only make use of V$ to compute the force, which suffices for our purpose (Syer & Tremaine 1996). 
We frequently refer to the Kalnajs function (Kalnajs 1971; Shu et al. 2000) which is explicitly defined by 

. 1 r[(m + 1/2 + ii.)/2]F[(m + 1/2 - \u)/2] (1/4 - iu/2) 
Nm{^) = K{u,m) = 2F[(m + 3/2 + i.)/2]F[(m + 3/2-i.)/2] = 2^ ' ^^^^ 

For m 7^ and setting ^ = vjvn, we write exp[i(m^ + v\aR}j\ = exp(imiy3). For m = 0, expression exp(i!/ln_R) should replace 
expression exp(im^) in treating axisymmetric perturbations (Syer & Tremaine 1996). We write perturbations in the form of 
exp[i(mS -I- i/lni?)] instead of exp(im(^) throughout our analysis unless otherwise stated. 

We note that from the scale- free condition (Lynden-Bell & Lemos 1999), there is no a priori reason that the scale-free 
index a (see equationUJ should be necessarily real. However, the imaginary part 9(q) will make such term as exp[i9(Q) ln_R] 
oscillating around 0. For example, the surface mass density E = i?" could be negative when Ss(q) 7^ 0. From the physical 
perspective, only a perturbation term may have such non-zero imaginary part of the scale-free index ex. 

Finally, we conclude that the gravitational potential of a scale-free disc is given by expression 11311 and the constraint for 
scale-free index /3 is /3 G (—1/2, 1/2) for cold discs (i.e., = 0) and /3 £ (—1/4, 1/2) for warm discs (i.e., k > 0). 



3 AN ISOPEDIC MAGNETIC FIELD ACROSS A SCALE-FREE THIN DISC 

The generation, amplification, and activities of magnetic field are widely believed to be caused by nonlinear MHD dynamo 
processes in an electrically conducting gas medium (e.g., Parker 1979; Moffatt 2000; Sofue et al. 1996; Beck et al. 1996). In 
contexts of star formation, numerical simulations (e.g., Nakano 1979; Lizano & Shu 1989; Basu & Mouschovias 1994) lend 
support to the notion that the mass-to-magnetic flux ratio A may be constant in both space and time. In Shu & Li (1997), 
A parameter was assumed to be a constant and two important theorems were derived based on this assumption. From the 
ideal nonlinear MHD equations, Lou & Wu (2005) demonstrated in a straightforward manner that a constant A is a natural 
consequence of the frozen-in condition for the magnetic flux. 

As in Section 2, the cylindrical coordinate system {R, 9, z) is adopted to describe the dynamics in an infinitely conducting 
razor-thin disc coincident with the z = Q plane. This razor-thin disc is threaded across by an 'isopedic' magnetic field that 
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exists for \z\ > in vacuum and is almost vertical to the disc plane at z = 0. As a natural consequence of the standard ideal 
nonlinear MHD equations (Lou & Wu 2005), we demonstrates the magnetic field to be isopedic (Shu & Li 1997; Shu et al. 
2000), with a spatially constant dimensionless ratio A of the mass per unit area T,{R, 6, t) to the magnetic flux per unit area 
Bz{R, 9, t) (perpendicular component of the magnetic field evaluated at 2 = 0), namely 

A = 27rG'^/^E/B^ = 27rG^/^A = constant , (16) 

where the factor makes A parameter dimensionless. 

Such kind of frozen-in isopedic magnetic field adds no extra complications to the hydrodynamic treatment of stationary 
perturbations in scale-free razor-thin discs for a global analysis. From the two MHD theorems established by Shu & Li (1997), 
we introduce the following transformation for an isopedically magnetized razor-thin scale-free disc 

$ ^ e$ , n ^ en , (17) 

where parameter e is the dilution factor for the disc self-gravity caused by the magnetic tension force, 

€=1- 1/A^ , (18) 
and G is the enhancement factor for the thermal gas pressure caused by the magnetic pressure force, 

e = l + (l + r,2)/(A^+77') . (19) 

In expression 11911 . parameter 77 is defined as the ratio of the horizontal gravity to the vertical gravity 27rGS (in absolute 
values) just above or below the disc plane, 

77 = |5lll/(2^GE) . (20) 

For an unmagnetized razor-thin scale-free gas disc, we simply have \ ^ 00 and set e = 1 and B = 1. If O is independent 
of R and 9 or if such dependence is nonessential and ignorable for simplicity, the replacement of 11 — > Oil is equivalent to 
H OH by definition Q . For an isolated scale-free disc sustained by its own self-gravity, the disc equilibrium would require 
|A| ^ 1 and thus e > 0. However, the existence of a massive background dark matter halo potential would allow jAj < 1 for a 
disc equilibrium and hence e < 0. This latter possibility leads to potentially interesting consequences. 

For a scale-free disc, it is possible to determine rj value more accurately. Considering an axisymmetric scale-free disc 
equilibrium, we take P G (—1/2, 1/2) and E — cR~^'^~^ with c being a constant coefficient. By the analysis of Qian (1992), 
we have the corresponding gravitational potential in the form of 

a>_ „^r.-2/3 7rr(l/2 - /3)r(/3) _ ^ct,-2I3y^(R\ 

* - "^^-^ r(i-/3)r(i/2 + /3) - "'^^ ^"^^^^ • 

The horizontal gravitational force g\\ tangential to the disc plane is given by 
511 = -V||$ = - (e^ A + $ = -2cGR'^"~^pY,{P)en . 

We therefore obtain 

Ifflll _ 2cGR-^^-^pY^{p) _ pYo{P) 
' ~ 27rGE 27rcGi?-2'3-i tt ' ^ ' 

leading to the following conclusions: r] € (0, -l-oo) when {3 G (—1/2, 1/2) and r) grows as /3 grows. Note that 77 = 1 when 
(3 = corresponds to a singular isothermal disc (SID; e.g., Shu et al. 2000; Lou 2002; Lou & Shen 2003; Shen & Lou 2004; 
Lou & Zou 2004, 2006). 

Effects of an isopedic magnetic field on a razor-thin scale-free disc are summarized below. The horizontal self-gravity 
force (?|| and the two-dimensional thermal gas pressure H of a scale-free disc will be effectively modified as 

9\\ ^ m\ ' n ^ en , (22) 

respectively. Or equivalently, we can express such MHD modifications as 

$ ^ e$ , H ~*<dH . (23) 



4 THE STATIONARY AXISYMMETRIC ROTATIONAL MHD EQUILIBRIUM 

We first specify the stationary axisymmetric rotational MHD equilibrium for a razor-thin scale-free disc with an isopedic 
magnetic field. This requires that every physical quantity does not depend on time t and azimuthal angle 9. In the following, 
we use subscript '0' to denote physical variables in the stationary axisymmetric equilibium. For a rotating disc, we also require 
vro = 0. Using these constraints in equations Q — HlOjl . we see that equations (O and l|10|l are satisfied automatically. The 
radial force balance @ appears in the form of 

vlo _ A{QHq + e<^Q + §o) 

R ^ dR ■ ^ ' 
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Naturally, condition 1241 represents that the centrifugal force, thermal pressure force, magnetic tension and pressure forces, 
the gravities from the disc and the background dark matter halo achieve a force equilibrium in the radial direction. The 
self-gravity potential "l>o and the enthalpy Ho are respectively replaced by e"I>o and 0i/o such that the effects of an isopedic 
magnetic field are included [see expression l|28^ ]. 

We first take Eo = coR~^^~^ and weo = boR~^ , where co and feo are two positive constants. Then the enthalpy becomes 

^0 = (f^^o" = m^kc?,^ R-'' , (25) 

and the corresponding self-gravity potential appears as 

$0 = ~GcoYoW)R-^'' . (26) 

Since the gravitational potential of the background dark matter halo is assumed to be axisymmetric {P = const) and 
unperturbed ($ = $o) in the presence of disc perturbations, we define the ratio of the background dark matter halo potential 
to the disc self-gravity potential as 

/ = l>/$o = f o/*o . (27) 

Ifere, / is a measure for the dark matter halo effect. The larger the / value, the more massive a dark matter halo is. 
The MHD radial force balance as described by equation H24|l then becomes 



boR = ^ = -2/3 



Ii±Mefecf - (e + /)Gc„y„(/3)' 
Zp 



i?-^"-' . (28) 



The background rotational MHD equilibrium simply requires the following relationship 

bl = 2l3GcoYoW){e + /) - (1 + 4/3)efccf . (29) 
According to definition (O, the sound speed Vso in the axisymmetric background disc is explicitly given by 

..^0 = fcFE-^ = il±Mfcef /t^+^^J?-" , (30) 

varying with R. A parameter w is defined as the square of the ratio of the sound speed Vso to the background rotation speed 
V0O, which is the inverse square of the rotational Mach number 

^-.|o-(f + 2/3) 6g • ^'^^ 
Physically, w is a measure for the disc temperature. A larger w value means that the disc temperature is higher, corresponding 
to a more random thermal motion in reference to the regular disc rotation. 

As a normalization, we set the radial self-gravity force in the axisymmetric disc be —1 at 7? = 1. Since the radial 
self-gravity force in the disc plane g\\ = —2f3GcoYo{l3)R~^'^~^eft, we then require 

2l3GcoYoW) = 1 . (32) 
Finally, the axisymmetric MHD disc equilibrium equation 12911 can be simply written as 

bl = . + f-bUl + 2m^ or bl^^^^il±IL_. (33) 

By the very definitions, both 69 and w should be positive. For a physical background dark matter halo potential, it is clear that 
/ ^ for an attractive background halo potential. We now demonstrate the physical requirement f + e ^ 0. When the disc 
is gravitationally isolated (/ = 0), we should have e 0, equivalent to A 55 1 (see eauation ll8|l . The physical interpretation is 
that the isopedic magnetic field should not be too strong to overtake the self-gravity in an isolated scale-free disc, otherwise 
the disc would be torn apart by the magnetic tension force. In the presence of a background dark matter halo potential, the 
tension force of the isopedic magnetic field is allowed to be stronger than the disc self-gravity to the limit of / + e > 0. The 
similar kind of effect can be seen for a composite disc system of two gravitationally coupled discs with an isopedic magnetic 
field, in which one is a stellar disc and the other is a magnetized gas disc. The stellar disc acts as the background halo 
potential, whose gravity counteracts the isopedic magnetic tension force in the gas disc (Lou & Wu 2005). 



5 DISPERSION RELATIONS FOR MHD DISC PERTURBATIONS 

We now introduce MHD perturbations in the axisymmetric equilibrium of a scale-free disc with subscript '1' revealing the 
association with the perturbation of a relevant physical variable. More specifically, we write 

E = Eo-HEi(i?,e»,t) , VR^VRi{R,e,t) , vg =vgo+veiiR,e,t) , $ = $0 + $i(ii,f,t) , $ = $0 , H = Ho + H^{R,e,t) , (34) 

respectively. Meanwhile, the magnetic flux should be also perturbed following the surface density perturbation. By the frozen- 
in condition, the ratio of the surface mass density to the magnetic flux remains unchanged at all times. That is, A remains 
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constant and thus e should be unperturbed (see equation 1181 . For disc perturbations, changes in ry would lead to variations 
in Q; following Shu et al. (2000), we presume that such induced variations of O are nonessential and are thus ignored. 

Substituting expression 13411 into the three basic MHD equations Q — IIOI I. we readily come to the linearized MHD 
perturbation equations in the forms of 

^ + l±iREovni) + 4s{^ovoi + ^iveo) = , (35) 



2veoV0i _ a(e//i + e$i) 



dt 



R 



R 



dvgi dvgo , Veo dvgi VgoVRi 



dR 



(36) 
(37) 



dt dR R de R Rde 

As both $1 and Hi can be expressed in terms of Ei together with the background equilibrium condition, we first take 

Ei(i?,6i,t) = ciR'^'^^~'^exp[i{ujt + me + u\nR)] = ciR~^'^' exp[i{LJt + m6)] , (38) 

where ci specifies a small amplitude coefficient and 

f3' = f3i- iu/2 . (39) 

In expression (I38II , exp(im0) and exp(i!/ln_R) represent the azimuthal and radial perturbations, respectively. As there is no 
a priori reason to force the scale-free index of the surface density perturbation to be the same as that of the background 
equilibrium, we allow /3i and (5 to be different, representing constant scale-free indices of the perturbation and background 
surface mass densities respectively. 

We then take both radial and azimuthal bulk flow velocity perturbations vri and wgi in the form of 

A{R) eyip[i{Lut + me)] , (40) 

where A{R) is a function of R to be specified by equations 1531 and 1541 for vm and vei, respectively. 
With equation 1381 and the Poisson integral, we obtain (Qian 1992) 

$1 = ~GciY^{f3')R-'"^' exp[i{Ljt + m9)] , (41) 

where /3' and Ym{u) are defined by equations 1391 and 1141 1. Properties of Ym{u) are summarized in Appendix A. 
For Hi being a perturbation of enthalpy H, we readily obtain 



Hi = fcrSo^^E 



,(1+4/3) -1/(1+2/3) D-2,3' r-/ , , n\] blwCi 2/3' , , /jm 

1 = k- -^Cq 'ciR ' eyip[i{Ljt + mO)\ = — R ' eyip[i{Ljt mO)\ 

(1 -f- 2/3) 



Co 



As /3i is only constrained by the surface mass distribution, we require f5i G (—1/2, 1/2) for either warm or cold discs. 
Here, the angular speed f2o and the epicyclic frequency «o of the axisymmetric background are given explicitly by 



— voo/R = boR 



-0-1 



2 
Kq 



RdQl/dR + ml = 2(1 -/3)ng =265(1-/3)7? 



-2/3-2 



After substituting expressions 13811 and 1401 into equations 1351 1 — 1371 . we readily obtain 

t(uj + mno)T,i + -^^(R^oVRi) H = , 

Jri OH rC 



i(uj -\- mQ,o)vRi — 2r2oi'ei 



i{Lo + mQ.o)vei + t^vri 
2i Iq 



djOHi + e^i) 
dR 

im{QHi + e$i) 
R 



Recombining equations 1451 and I46II . we obtain 



and 



vri = 



vgi 



[{uj + mflo) 



(eiJi + e$i) 



(ei7i-fe$i) 



Substituting expressions 13811 . 14H . 1421 . 1471 and 1481 into equation 1441 . we immediately arrive at 



{uj + mQ.o)R-^^' + [wQhl - ecoGYmiP')] ^ 



R 



-2/3 



[wQbl - ecoGYrnip')] 



dR \ [{u + mno)^ 
mR-^^-^ 



R 



-2/3' 



, ^ d 



R'^^' = 



(42) 

(43) 
(44) 
(45) 
(46) 

(47) 

(48) 



(49) 
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where Q.o and are defined in equation I43II . As the angular frequency is a constant, we derive a simpler form of equation 
1491 1 by writing uj = ujR^'^^ , namely 

{o + mbo _ mho K + 2{I3 + 13'- 2(3'^)] + tl^ [m^ - 2/3' (2/3' - 1)] 



wQhl - ecoGYr^m {uj + m6o)2 - 2(1 - (3)hl 

4(1 + /3) [mho - P'ico + mho)] ii(tl' + mho) 
^ [(^i + mbo)2-2(l-/3)fo2]2 ■ 

Equation 15011 can be viewed as the local dispersion relation for time-dependent MHD perturbations in a scale-free razor-thin 
disc because of the explicit 7?— dependence of cj. Apparently, it is a complex quintic equation in terms of uj. According to 
the Abel Impossibility Theorem^, we cannot solve for lj analytically from equation l|5()|l generally. Nevertheless, this does not 
prevent straightforward numerical solutions. 

The classical Wentzel-Kramers-Brillouin- Jeffreys (WKBJ) relationship in the so-called tight-winding regime can be readily 
recovered with certain simplifications and approximations. In the parameter regime of Qw <^ 1, v ^ m and cu 0, dispersion 
relation ^bOl naturally reduces to 

((i + mhof = 2(1 - l3)bl + wehlu"^ - 2-K\u\ecoG . (51) 

By defining the radial wavenumber k as v/R, the WKBJ dispersion relation I51II can be cast into the form of 

{u + mQ.of = Ko + k^Qvlo - 2TieT.oG\k\ , (52) 

(see Lou & Fan 1998a for fast MHD density waves) which is essentially the same as equation (39) of Shu et al. (2000). This 
consistent correspondence to the classical WKBJ dispersion relation is a necessary check. Equation 1)50^ can be viewed as 
a more comprehensive local description for MHD perturbations in a scale-free disc, in comparison to the standard WKBJ 
approximation. Detailed derivation procedures and formulae can be found in both Appendices B and A. 

By equation 11811 . e is a function of A and by equations 1191 and 1211 1. is a function of both A and 13. From equation 1321 
that normalizes the tangential gravity force in the disc plane, we determine cq by a specified (3 value. According to equation 
l|33|l . ho parameter can be expressed in terms of /, w, (3 and A. We thus formally write out the solution of dispersion relation 
H5()|l as = T{m, V, f3i; f3, f, w, A), yielding u) — R~'^~^J^{m, v, /3i; f3, f, w, A), where parameters m (an integer for 
the azimuthal wavenumber) , f (for the radial wavenumber) and /3i (the scale-free index for perturbations) characterize MHD 
perturbations, while parameters (3 (the scale- free index for the background), / (the indicator for the background dark matter 
halo), w (the indicator for the disc temperature) and A (the indicator for the magnetic field strength) specify the axisymmetric 
background rotational MHD equilibrium. As dispersion relation H5()|l is a complex polynomial, MHD perturbations in a disc 
are stable only if or .7-" is real. 

In order to derive corresponding forms of vm and vei from equation H40|l . we directly substitute expressions 1)38^ . I|411 
and 1421 into equations 1471 and 1481 to explicitly obtain 

VRi = — [wQho - ecoGFm(/3 . -7- ^TT^T-^ exp i(Ljt + me)] , (53) 

Co [(mho + ui)'' — 2(1 — f3)hg\ 

Because uj in expressions 153|l and H54|l contains the i?— dependence of a) = ujR'^'^^, both velocity perturbation components 
VRi and vgi no longer take a simple power-law form in R besides the wave factor exp[i(Ljt + m9 + vlnR)]. 

A special subcase of equation l|50^ is the axisymmetric case with m = 0, from which we immediately have 

1 2(3' (2/3' + 2/3+ l)uj^ - 4^13' (2/3' - 1) (1 - /3)&g ' 



uebg - ecoGFo(/3') ' [<i2 - 2(1 - /3)b2]2 < °- ^^^^ 

Equation 1551 can be solved analytically and obviously cj = is one of the solutions. In discussing the axisymmetric case, we 
therefore adopt uj ~* instead of uj — 0, for the stationary condition in order to avoid the singularity caused by lj = 0. 

Having admitted our limited capability to study time-dependent local dispersion relation 1501 analytically, we now explore 
the global stationary perturbation configurations, which is simpler yet important (Shu et al. 2000; Lou 2002; Shen & Lou 
2004a; Shen, Liu & Lou 2005; Lou & Bai 2006). Stationary MHD perturbation configurations or zero- frequency neutral MHD 
perturbation modes correspond to lj = in dispersion relation H50^ . For clarity, we emphatically distinguish the axisymmetric 
and non-axisymmetric cases at this point. 

For non-ajcisymmetric cases with m 7^ 0, we immediately have 

hi _m' + 2(/3 + 13' - 2/3'^) 



wQhl-ecoGYrr,((3') m-^- 2(1 -13) 



(56) 



^ A general polynomial equation of order higher than the fourth degree cannot be solved algebraically in terms of finite additions, 
multiplications and root extractions (e.g., Hungerfort 1997). 
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With normalization I32II . we readily cast equation 1561 into the following form of 

- 2(1 - 13) 



m2 + 2(/3 + /3' - 2/3'2)J ' ^^'^^ 

Because of the symmetry in m for Ym{(3') = Y-m.{P') as shown in Appendix A, equation I56II is valid for both m > and 
m < where m is an integer. From now on, we focus on m ^ without loss of generality. 

The axisymmetric case of m = is somewhat special. We easily see the axisymmetric dispersion relation 1)55^ is satisfied 
for uj = 0, because zero is a valid solution for uj when m = 0. Alternatively, we let ui approach zero {ui —> 0) but do not vanish 
(oj 7^ 0). We adopt this limiting procedure to deal with the axisymmetric case (Lou & Shen 2003; Shen & Lou 2004b; Lou 
& Zou 2004, 2006; Lou & Wu 2005; Lou & Bai 2006). Using this axisymmetric limiting procedure, we determine the global 
stationary axisymmetric MHD dispersion relation in the form of 

bl _ I3'{2f5' - 1) 



[™efo2 _ ecoGFo(/3')] (1-/3) 
Using normalization 13211 . dispersion relation 1581 can be rewritten as 



2f3Yo{l3) 



(58) 



(59) 



f3'{l - 2/3') 

Equations 15711 and I59II appear similar to the nonaxisymmetric self-consistent relation (39) and the axisymmetric self- 
consistent relation (57) of Syer & Tremaine (1996) respectively except for the two extra parameters e and Q representing the 
effect of an isopedic magnetic field and for a difi^erently defined /3' parameter. In the absence of an isopedic magnetic field 
with e — 1 and O = 1 and with the perturbation scale-free index /3i being the same as the background scale-free index f3, 
equations 1571 and 1591 1 respectively reduce to equation (39) and (57) of Syer & Tremaine (1996) precisely. 

In the analysis of Syer & Tremaine (1996), the global aligned and logarithmic spiral perturbation configurations are 
integrated together by the parameter ly, i.e., u = for aligned cases and v ^ for logarithmic spiral cases. However in their 
work, the axisymmetric and nonaxisymmetric cases are handled in different approaches. Here as a unification in our analysis, 
equation l|50|l is equally valid for axisymmetric (m = 0) and nonaxisymmetric (m 7^ 0) perturbations, as well as for aligned 
{u — 0) and logarithmic spiral perturbations (z^ 7^ 0). Therefore, equation 1501 becomes a very general and rigorous dispersion 
relation for global stationary MHD perturbation configurations and for axisymmetric stability of a MHD scale-free disc. 

Under the stationarity condition, the two velocity components vm and vei can be written as 

c, [wObl - ecoGYmiP')] 2im(l-/3') ^.2p'+0 
^ ^0 K -2(1-/3)] '^ eM^me), (60) 

^ K- 2(1-/3)] ''exp(rme). (61) 

When m = 0, we have vr\ — according to expression 1601 indicating no radial velocity perturbation. 

By convention, the real parts of expressions I6OII and 16H represent the physical solutions for vri and voi- In addition 
to the wave perturbation factor exp[i{u!t + mO -\- fin J?)], amplitudes of vr\ and vg\ carry the power-law of R in the form of 
^-2/3i+/3^ For /3i = /3, perturbation velocity components vri and vg\ share the same power law of R as R~^ , similar to the 
background power law in R such as expression Ulljl . 

In general, it would be natural to take /3i = (3, which can be explained by a direct proportionality between a MHD 
perturbation and the background equilibrium. However, for logarithmic spiral perturbations, /3i = 1/4 can insure a constant 
radial flux of angular momentum (Goldreich & Tremaine 1979; Syer & Tremaine 1996; Fan & Lou 1997; Shen & Lou 2004a). 
For aligned cases, there is no angular momentum flux inwards or outwards. For logarithmic spiral cases, the total radial 
angular momentum flux (including the gravity torque Ag, advective transport Aa and magnetic torque Ab) is given by the 
sum of Atotai = Ag + Aa + Ab oc T?^"**^^. Therefore, only for /3i — 1/4 can the radial angular momentum flux be independent 
of 7?. Detailed procedures can be found in Appendix G. 

Although dispersion relation 15011 is derived for MHD perturbations without the WKBJ approximation, it should be 
properly viewed as a local dispersion relation for lu 0, otherwise lu cannot remain constant but varies with R. As such, 
equation H50|l is complex in general. Only under two special circumstances can dispersion relation H5()|l become real. One is 
when u — 0, while the other is when lu — and /3i = 1/4. The flrst case of 1/ — means no radial wave oscillations and thus 
no radial flux of angular momentum, while the second case of lj = and /3i = 1/4 assures a constant radial flux of angular 
momentum (see Appendix G). Furthermore, dispersion relation 15011 should be regarded as a global stationary dispersion 
relation when a; = 0. It might be suggestive of a link between whether the dispersion relation 15011 is real or complex and 
whether the radial angular momentum flux is constant or not, as no angular momentum flux {1/ — 0) can be seen as a special 
case of a constant radial flux of angular momentum in R. It seems that a non-constant radial flux of angular momentum 
makes dispersion relation l|50|l complex. It is a common choice to assume a constant angular momentum flux in R (Goldreich 
& Tremaine 1979; Lemos et al. 1991; Shu et al. 2000; Lou & Shen 2003; Shen & Lou 2004a; Lou & Wu 2005). 

By imposing the stationary condition of 0, dispersion relation 1501 represents a constraint on the background 
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parameters (/, w, /3, A) and the perturbation parameters (m, u, simultaneously. While a real form of equation 15UII gives 
one constraint, a complex form would then give two. The additional constraint is due to the imaginary part of equation (Syer 
& Tremaine 1996), which gives a more strict condition on possible parameters (see Figs. 2a and 2b below). 

In the following analysis, we mainly focus on global stationary MHD perturbation configurations. Since the O parameter 
appears always together with w in the dispersion relations, we introduce a new variable combination w = Ow for simplicity. 
A natural interpretation for this can be related to the very definition of w for Vso/vgQ. In other words, w can be defined as 
{VQvbo)'^ / veo with ^/QvsO representing the magnetosonic speed (Lou & Wu 2005) and seen as a measure of disc temperature 
and magnetic pressure together. We emphasize here that by equation l|19|l . the value of Q parameter can be very large as 
— > and A — > 0. In disc galaxies, effects of magnetic pressure and tension can be limited in a certain manner under specific 
situations. From now on, we shall deal with w — Qw throughout the remaining analysis. 

Obviously, we should require ^ on the ground of physics. As a repulsive background dark matter halo potential 
is unphysical, we also demand / ^ 0. According to equation 1331 1. we need to impose / + e ^ to maintain a background 
rotational MHD equilibrium of axisymmetry. In constructing global stationary MHD perturbation configurations, we keep in 
mind these three constraints, namely, w ^ 0, f ^ and / + e J? 0. 

For the perturbation scale-free index we have /3i G (—1/2, 1/2), while for the background scale-free index /?, we 
have /3 G ( — 1/2, 1/2) for cold discs and f3 € (—1/4, 1/2) for warm discs. But for the convenience of analysis, we still use 
P G ( — 1/2, 1/2) for general situations hereafter. One should keep in mind that when (3 G (—1/2, —1/4], the variable w can 
only take on zero value. 



6 GLOBAL STATIONARY ALIGNED CASES WITH = 

For aligned global MHD perturbation configurations, the isodensity contours of perturbations are aligned in azimuth at various 
radii. Physically, such aligned cases should be interpreted as purely azimuthal propagation of MHD density waves (Lou 2002). 
Without the constraint on the angular momentum flux transport and as an example of illustration, we may take f3i = (3 
indicating that MHD perturbations carry the same scale-free index as the background one. 



6.1 Axisymmetric Disturbances with m — 

For axisymmetric configurations with m = 0, expressions 13811 . 1601 and 16111 give the following results 

El = ciR ' , Vri = , vgi ^ ^ — - — ; -R , (62) 

2co&o 

which are trivial in the sense of a perturbation rescaling for the axisymmetric background rotational MHD equilibrium (Shu 
et al. 2000; Lou 2002; Lou & Shen 2003; Lou & Zou 2004, 2006; Lou & Wu 2005). 



6.2 Nonaxisymmetric Disturbances with m > 

When m > 0, equation 157|l can be readily rearranged into the form of 



n? - 2(1-/3) 



ie + f) 



ra^ - 2(1 - f3) 



m2+ 4/3(1-/3) 



(63) 



[1 (l-f 2/3)w] 

where i;; = Ow and 69 is determined by equation 1331 . For the range of /3 G (—1/2, 1/2), numerical computations show that 
2/3Yo(/3) remains always positive after taking limit at certain points, such as /3 = 0. Furthermore, with the asymptotic form 
of Ym{fi) —> 2-71 /m as m — > 00, it is guaranteed that Ym{f3) > for m 1 based on carefully designed numerical tests. 
For simplicity of expressions, we introduce a handy parameter 

[m^ - 2(1 - (3)] 



Ca 



(64) 



[m2-h 4/3(1-/3)] ■ 

When the scale-free index /3, the azimuthal wavenumber m and the gravity dilution factor e are specified, equation H63^ 
establishes a unique relation between the inverse square of the magnetosonic Mach number (denoted by w parameter) and 
the background gravitational potential ratio (represented by / parameter) that sustains a global nonaxisymmetric stationary 
MHD perturbation pattern. 

A more straightforward relation between / and w can be simply expressed as 
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Figure 1. (a) The physically allowed combinations of w (a measure for the disc temperature and magnetic pressure together) and / 
(a measure for the background dark matter halo potential) for global stationary MHD perturbation configurations with the scale-free 
index /3 = —0.22 and e = 0.5 (i.e., the disc self-gravity overtakes the magnetic tension force); the two solid curves separately stand for 
m = 1 and m = 2 aligned cases, while the dotted curve stands for m = 1 unaligned case and the curves of m ^ 2 unaligned cases are 
unphysical. Note that different points on the dotted curve have different v values and the bifurcation point is at (461.2, 0.2764) in the 
(i2), /) diagram, (b) The same as in panel (a), except for e = —0.5 (i.e., the magnetic tension force overtakes the disc self-gravity) and 
thus the requirement / > 0.5. The unaligned logarithmic spiral cases all become unphysical. 



where Ci and C2 are defined explicitly by 

= — wm — ' ^ ^ 

When e = corresponding to the situation of the magnetic tension exactly cancelled by the horizontal disc self-gravity, 
equation 15(il l and background equilibrium relation 13311 then give 

W = Ca, f = hl[l + {l + 2l3)Ca\ ■ (68) 
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In order to ensure Ca 0, we need m 2 or /3 G (—1/2, —0.2071). Thus for e = 0, aligned MHD perturbation configurations 
of m 2 can exist, while the m = 1 mode can only exist when /3 G (—1/2, —0.2071). As the horizontal disc self-gravity is 
cancelled by the magnetic tension force, only the dark matter halo potential resists the centrifugal force, the thermal pressure 
force and the magnetic pressure force to maintain the disc equilibrium. Note that when e = 0, we have a fixed w and bg is 
determined by / as compared to the expression of bg using w and / in equation H'A'A^ . This situation is strikingly similar to 
those of low-mass discs (/ — > oo) referred to by Syer & Tremaine (1996). Such a similarity can be easily understood because 
when e = 0, the horizontal disc self-gravity is completely cancelled by the magnetic tension force and the effective / (defined 
as the ratio of the force from the background halo to the force from the disc) apparently approaches infinity. In a certain 
sense, e = represents a critical point separating two qualitatively different cases. For < e ^ 1, the magnetic tension force 
is weaker than the horizontal self-gravity force of the disc, while for e < 0, the magnetic tension force becomes stronger than 
the horizontal self-gravity of the disc. One would expect systems dominated by magnetic Lorentz force or gravity to show 
different physical properties. 

When < e 5C 1, the entire analysis here parallels the corresponding part of aligned discs in Syer & Tremaine (1996), 
because we can simply take f/e and Qw to replace their / and the w variables, respectively. In other words, the two factors 
1/e and Q can be regarded as linear rescaling factors of / and w variables when the magnetic tension force is weaker than 
the horizontal disc self-gravity force. 

When e < 0, the magnetic tension force becomes stronger than the horizontal self-gravity of the disc. From the background 
rotational MHD equilibrium condition 13311 . the gravitational potential of a dark matter halo must exist and be strong enough 
in order to create an MHD disc equilibrium in the first place. As already noted earlier, physical constraints on the disc system 
are lo ^ and / ^ max{0, — e}. With these requirements, we can infer the physically allowed part of a hyperbolic curve of 
equation 1)6 5 |l in the {w, f) diagram and determine the relevant range of w. 

For the aligned cases, we reach the following conclusions drawn from the analysis detailed in Appendix A. For m 2 
MHD perturbation modes (bar-like), we have w G (Ca , Ca + C1/C2] when < e ^ 1 (see Figs, la, 2a, 2b) and w £ [0 , Ca) 
when e < (see Figs, lb, 2c). For m = 1 MHD perturbation modes with < e ^ 1, we have w G (Ca , +00) when 
/? G (-1/2 , -0.2071) (see Fig. la), w G [C1/C2 +Ca , +00) when /3 G (-0.2071 , 0) (see Fig. 2a), and m G [0 , C1/C2 +Ca] 
when /3 G (0 , 1/2) (see Fig. 2b). For m = 1 MHD perturbation modes with e < 0, only when /3 G (-1/2 , -0.2071) and 
w G [0 , Ca) can global stationary MHD perturbation patterns possibly exist (see Fig. lb). 

We now discuss several special cases below. 

The case of /3 = corresponds to a flat rotation curve of constant vgo and a singular isothermal disc (SID) under scale- free 
conditions (Shu et al. 2000; Lou 2002; Lou & Shen 2003; Shen & Lou 2004a; Lou & Zou 2004, 2006; Lou & Wu 2005; Lou & 
Bai 2006). Using the asymptotic expression r(a;) ~ in the limit of a; ^ 0, we write equation l)65|l as 

^ e(m — l)(mw — m — 2) e(m — 1) 
(m'^w — m? + 2) m 

where /? = has been set. One can readily see that / = for an arbitrary w when m = 1 and l3 — 0. In other words, an 
isopedically magnetized isolated disc with a flat rotation curve supports an m = 1 MHD perturbation mode for an arbitrary 
Mach number (Syer & Tremaine 1996; Shu et al. 2000; Lou 2002; Lou & Shen 2003; Lou & Zou 2004, 2006; Lou & Wu 2005). 
For m ^ 2 MHD perturbation modes, we have w G (1 — 2/m^ , 1 + 2/m] when < e ^ 1 and w G [ , 1 — 2/m?) when e < 0, 
respectively. 

Another singular case appears when /3 = (1 — \f2)/2 — —0.2071 and m = 1. According to equation H56|l . we have 60 = 
corresponding to a nonrotating disc. In other words, for a disc with a scale-free index /3 = —0.2071, it has to be nonrotational 
in order to support a stationary m = 1 MHD perturbation mode. 



2(m + 1) 

m?w — (m? — 2) 



- 1 



(69) 



7 GLOBAL STATIONARY UNALIGNED LOGARITHMIC SPIRALS WITH v^Q 

For unaligned logarithmic spiral cases, the isodensity contours of MHD perturbations entail a systematic phase shift in azimuth 
as R increases, which gives a logarithmic spiral pattern. Unaligned stationary logarithmic spirals involve both azimuthal and 
radial propagations of MHD density waves. For a constant radial flux of angular momentum, we require /3i = 5R(/3') = 1/4. 
To be general, we take /3i and (3 to be two different parameters. 



7.1 Nonaxisymmetric Disturbances with m 7^ 

For m 7^ 0, the corresponding relation 15711 can be written as 
eB^bl{w~C) , 



(70) 



14 Y. Wu & Y.-Q. Lou 



where 

^ = ^) = ^) ' ^ ^ [^^r2(^f^"-V)] = ■ ('^^ 

Here, parameters bg, w and e are all real, while functions B and C may be complex in general. We then have 

e?ft{B) = &o[u' - 3^(C)] , e9(B) = -bl^iC) . (72) 

We can also write the real and imaginary parts of C, 5R(C) and 9(C), explicitly as 

^f^. ^ W - 2(1 - /3)](m^ + 2/3 + 2/3i ^ 4/3? + t/^) 

^ (m2 + 2/3 + 2/3i-4/32+jy2)2 + j,2(i_4^^)2 ' ^'•'i 

c^rc) ^ Kl-4/3i)[m^- 2(1-/3)] 

(m2 + 2/3 + 2/3i- 4/32+ 1.2)2 + j,2(i_ )2 • W^tj 

Based on the analysis of Appendix A and equations 17311 and I74I I. we notice the following relations that C(m, —v, /3i; fH) — 
C*(m, u, /3i; /3), B(m, -i/, /3i; /3) = B*(m, i^, /3i; /3), B(m, i^, /3i; /3) = B(-m, i/, /3i; /3) and C(m, i/, /3i; /3) = 
C(— m, i., /3i; /3), where the asterisk * indicates the complex conjugate operation. Therefore, equation 17211 remains valid for 
±1/ and ±m. This conclusion can be seen as a manifestation of the anti-spiral theorem (Lynden-Bell & Ostriker 1967), which 
states that trailing spiral and leading spiral share the same solution forms under stationary and time-reversible conditions. 
Therefore, there is no loss of generality to consider m ^ and 1/ > 0. 

When P G (-1/2, 1/2) and m ^ 1, it is clear that - 2(1 - /3) / 0. According to equation ifTIjl . we thus have 9(C) = 
only if u — or f3i = 1/4. Obviously, 1/ — corresponds to the aligned case as has been discussed already. The requirement 
/3i = 1/4 implies a constant radial flux of angular momentum associated with logarithmic spiral MHD perturbations. 

When e = 0, we need /3i = 1/4 and m ^ 2 to satisfy equation H72|l . 



7.1.1 Discs with [3i = 1/4 

When /3i — 1/4, equation 1701 becomes real and appears as 

r{m/2 + l/4 + iiy/2) 



2/3^0 (/3) 



r(m/2-H3/4 + ii//2) 



[1 + {1 + 2P)w] 



- 2(1 - /3) 



(m2 2/3 1/4 + 1/2) 



Do 



Di 



which can be further transformed into 

/ = 

where the three coefficients Do, Di and D2 are explicitly deflned by 
[m^ - 2(1 - /3)] 



or (f + eD2)iw ~ Do) = eDi , 



(75) 



(76) 



Do 
Di 
D2 



(m2-f 2/3 + 1/4-f 1/2) 



2/3^0 (/3) 
, 7r(l - 



r(m/2 + l/4 + ii//2) 



r(m/2 + 3/4 + ii//2) 



[2m2(l + /3)+ 4/32 -7/4 + 1/2] 



(m2 + 2/3 + l/4 + !/2) 



2/3) 



2/3Fo(/3) 



r(m/2 + l/4 + ii//2) 



r(m/2 + 3/4 + ii//2) 



A relationship is thus established by equation 1761 between w and /. By specifying parameters m, u and /3, this relation 
gives rise to a hyperbolic curve in the {w, f) diagram. From Appendix B, we thus draw the following conclusions. For m ^ 2 
MHD perturbation modes, we have w G (Do , Do + D1/D2] when < e ^ 1 (see Fig. 3a) and w £ [0 , Do) when e < (see 
Fig. 3b); for m = 1 MHD perturbation modes, we have w £ [0 , Do + D1/D2] when < e ^ 1 (see Fig. 3a) and these modes 
do not exist when e < 0. 

When e = 0, we have a fixed w as to = Do and thus m = 1 MHD perturbation mode cannot exist for a negative Do. A 
relationship between / and feg can be established in the form of 



/ = 6g [1 + (1 + 2/3)Do] 



(77) 



which is similar to equation I68II for the aligned cases. 
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Figure 3. (a) Given the baclcground scale-free index /3 = 1/4 and e = 0.5 (i.e., the horizontal disc self-gravity overtakes the magnetic 
tension force), the physically allowed combinations of id (a measure of the disc temperature and magnetic pressure together) and / (a 
measure of the background dark matter halo potential) for m = 1 and m = 2 global stationary MHD perturbation configurations of 
the perturbation scale-free index /3i = 1/4 disc with f = (aligned) and u = 5 (unaligned), respectively. Two solid curves stand for 
m = 2 MHD perturbation mode and two dotted curves stand for m = 1 MHD perturbation mode, (b) The same as panel (a), except for 
e = —0.5 (i.e., the magnetic tension force overtakes the horizontal disc self-gravity). Note that a m = 1 MHD perturbation mode does 
not exist for e < 0. 



7.1.2 Stationary Logarithmic Spiral Patterns in Discs with /3i 7^ 1/4 



When /3i 7^ 1/4, we have Si(C) 7^ 0. Equation H7()^ is then complex, giving rise to two real equations. From equation H72^ . we 
solve for w and / to obtain 



3ff(C) 



9(B) 
1 + (1 + 2/3) (5R(C) 



3f?(g)9(C) 
9(5) 



+ 1 



(78) 
(79) 
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From extensive numerical computations, we gain the following useful information. First, 5R(-B) > and decreases with 
increasing ly. Secondly, 9(B) > for -1/2 < /3i < 1/4 and Q{B) < for 1/4 < jSi < 1/2, respectively. 

We emphasize here that the w and / solutions (1781 and (I79II give a specific point in the {w, /) diagram, which differs from 
the aligned case with i^ = 0or/3i = l/4 discs whose real dispersion relation gives a hyperbolic curve in the (w,/) diagram. 
Under the conditions u and /3i 7^ 1/4 (non-constant radial flux of angular momentum, see Appendix G), dispersion 
relation l|70|l is generally complex, corresponding to two constraints (i.e., real and imaginary parts of equation I70|l on the 
combination of w and / and reducing a solution curve to a solution point in the {w, /) diagram (see Figs, la, 2a, 2b). 

Finally, we reach the following conclusions numerically: ^i; < when m ^ 2; for m = 1 MHD perturbation modes, w 
remains always positive, while / > only for e > 0. Therefore, global stationary spiral MHD perturbation configurations with 
/9i 7^ 1/4 can only exist for m — 1 MHD perturbation modes and when e > (see Figs, la, 2a, 2b). 



7.1.3 Bifurcation Points from Aligned Cases 

The sequence of stationary logarithmic spirals with u G (0, +00) bifurcates from the aligned discs with = at the 
critical point where 0. When v becomes small, Ss(-B) and Q'(C) may be replaced by —v{dB/dp')/2 and —v{dC/dp')/2 
respectively, where the two first derivatives are evaluated in the limit of ^ (/3' = /9i). 

Therefore at the bifurcation point (indicated by a subscript j-|j£ of a physical variable), equations 17211 become 

eB^bl[w-C), e^ = -hl^, (80) 

where all quantities are evaluated at i/ = {(3' — I3\). Eliminating 6q between the above two expressions, we obtain 



™bif = 



( dC\ f dlnB \-^ 



where the derivatives are evaluated at = or = 
We then obtain 

2(l-4/3i)[m^ + 2(/3-l)] [m^ +2(^-1)] 

bif [(m2 + 2/3 + 2/3i - 401)^^^(0,)] ^ (m^ + 2/3 + 2/3i - 4/3?) ' > 

where 



^2 2J ^ \ 2 2 

and ^(z) is the digamma function (e.g., Shen & Lou 2004). Once w is determined, the expression of / follows from equation 
1651 1. namely 



n . ^o, , (m^ + 2/3 + 2/3i-4/3?)[m^(l + /3)-(l-/3i) + 2(/3^-/3i^)] ,^, 
+ ^f^' + (i_4/30K+2(/3-l)l ^'"(^^^ 



(82) 



^bif \2(3YoW) 

Together, this pair of w^^jf (eQuation l81ll and f-^-^^ feauation l82ll gives a specific point in the (w), /) diagram which signifies the 
bifurcation place between the aligned and unaligned sequences (see Figs, la, 2a, 2b) . Numerical explorations have demonstrated 
that when m ^ 2, w^^jj- < 0, while for m — 1, both w^-^^ and f^-^^ are positive. Therefore only for m = 1 MHD perturbation 
modes, the bifurcation point is physically meaningful. 



7.2 Marginal Instability of Axisymmetric Disturbances (m = 0) 

When m = 0, dispersion relation 15911 can be written as 



= bo 



2/3Yo{f3) /3i(l-2/3i) + //V2-i(l-4/3i)///2 

We adopt Pi — 1/4 to make equation I83II real, namely 



w + 



(1-/3) 



(1 + 4z/2) 



where A/'o(i') is the Kalnajs function defined by equation HSU . From equation I84I I. we then obtain w as 



(1 + 2/3) 



8Yo(/3)/3(l-/3)(l + //6) 



7r(l + 4Z/2) 



/3yo(/3)(l + .//e) 



(83) 



(84) 



(85) 



7r(l + 2/3) 

where is defined as 1/w = l/(Ow), which is the same as the introduced by Shu et al. (2000). When e < 0, we requne 
e + / > and thus 1 + //e < 0. Therefore when e < with 1 + f /e < 0, w remains always negative according to equation 
1851 1. When e = in equation 18411 . there is no global stationary configuration, as in the situation of e < 0. 
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Figure 4. The axisymmetric (m = 0) marginal instability curves of D-^ i^/w) (the reciprocal of the disc temperature and magnetic 
pressure together) versus u (the dimensionless radial wavenumber) with different ratio f/e equal to (dashed curves) and 0.2 (solid 
curves), respectively, for specified values of the perturbation scale-free index /3i = 1/4 and the background scale-free index (3 = 0. We 
take e > (i.e., the disc self-gravity overtakes the magnetic tension force). For e < (i.e., the magnetic tension force dominates the 
disc self-gravity), the instability does not exist. Note that both the Joans collapse and ring fragmentation unstable regions shrink with 
increasing f/e. 



When equation 18511 is displayed in the {ly, D^) diagram, the curves correspond to the marginal instability {lo = 0). Here 
u represents the radial wavenumber of MHD perturbations, and a smaller u signifies a larger perturbation scale. Also, is 
the square of the rotational Mach number, and a larger thus signifies a faster rotation. Using equation 15211 . we can analyze 
stability properties of the {u, D^) diagram in different parameter regimes separated by the marginal instability curves (Shu 
et al. 2000; Lou 2002; Lou & Shen 2003; Shen & Lou 2004b; Lou & Wu 2005; Lou & Zou 2006). Here uP > Q means stability 
and tj^ < means instability, respectively. According to equation H52|l . we see that when e <0, uP remains always positive so 
that the disc is completely linearly stable to axisymmetric perturbations at all wavelengths. As a consequence, the marginal 
stability curve in Figure 4 does not exist when e < 0; in other words, both Jeans collapse and ring fragmentation instabilities 
are suppressed (i.e., lower-left curves become negative and upper-right curves move to infinity). 

For e > 0, we discuss only a special case of /3 = as an example of illustration, which corresponds to the SID case (Shu 
et al. 2000; Lou 2002; Lou & Shu 2003; Lou & Wu 2005). A more general treatment (Shen & Lou 2004a, b) wiU be given in 
a forthcoming paper especially aiming at the axisymmetric instability. 

Knowing the asymptotic limiting expression I3Yq{I3) ^ tt as /3 0, we obtain l/w as 

i = n2 ^ (l + /A)--^oM ^86) 

w- AAo(/.)-8(l + /A)/(l + 4/.2) ■ 

From Appendix F, the solution curves of equation 1861 separate the {v, D^) diagram into three or two regimes. The lower- 
left corner and the upper-right corner are the unstable regions with lP < 0. The lower-left corner corresponding to MHD 
perturbations of large radial spatial scales represents the region of Jeans collapse instability modified by the disc rotation 
and the isopedic magnetic field, while the upper-right corner corresponding to perturbations of smaller radial spatial scales 
represents the region of the MHD ring fragmentation instability (Lou & Fan 1998a; Shu et al. 2000; Lou 2002). Sandwiched 
between these two unstable regimes is the stable region of oj^ > (Shen & Lou 2004a; see Fig. 4 for direct visual impressions). 

We can also demonstrate that as / increases or e decreases, both the MHD Jeans collapse and ring fragmentation 
instability regions tend to shrink. In other words, an enhanced background dark matter halo potential or a stronger isopedic 
magnetic field tend to squash instabilities both in small ( MHD ring fragmentation) and large (MHD Jeans collapse) spatial 
scales. In particular, when 1 -I- f/e > No{0) = 4.337, the MHD Jeans collapse instability region disappears. In contrast, the 
MHD ring fragmentation region always exists. For example, a cold disc (?j) = or +oo) will always subject to the ring 

fragmentation instability. Our numerical explorations are supportive of some conclusions of numerical simulation results of 
Miller (1978), namely, a certain level of velocity dispersion is necessarily required for the disc stability, and furthermore, a 
massive dark matter halo cannot stabilize a 'cold' disc galaxy completely. 

In short, we come to the conclusion that a massive background dark matter halo (i.e., a high / ratio) and a strong isopedic 
magnetic field (i.e., a small e) tend to make the magnetized gas disc more stable against axisymmetric MHD disturbances. 
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8 CONCLUSIONS, SUMMARY AND DISCUSSION 

Using the standard linearization procedure for MHD perturbations in a rotating disc, we have in a more general manner 
constructed global stationary MHD perturbation configurations (i.e., small deviations from the axisymmetric background) for 
both aligned and unaligned logarithmic spiral patterns in a razor-thin scale-free disc system embedded in an axisymmetric 
dark matter halo potential and with an isopedic magnetic field almost vertically threaded through the gas disc plane. The 
stability analysis of axisymmetric MHD perturbations is performed in a disc of flat rotation curve with /3 = (i.e., SID) as 
an example of illustration. The most important and interesting contribution of our paper is that we integrate several essential 
aspects, such as isopedic magnetic field and dark matter halo, into the disc dynamics and adopt an analytic approach for the 
stability analysis. A generalized relationship is achieved and some interesting properties directly derived from the relationship 
have been revealed. 

As global stationary MHD perturbation configurations may represent transitions between stabilities to instabilities and 
vice versa (Shu et al. 2000), we can use the obtained global MHD perturbation configurations to explore initial conditions for 
possible instabilities and the dynamical evolution of a razor-thin scale-free disc system. Thus, a more realistic magnetized disc 
system to model actual spiral galaxies may provide sensible physical information. Since the systematic galactic observations 
by Rubin et al. (1982) and Kent (1986, 1987, 1988), the massive dark matter halo becomes an indispensable component for 
models of disc galaxies, which are referred to as a partial disc system (e.g., Binney & Tremaine 1987; Syer & Tremaine 1996; 
Shu et al. 2000; Lou & Shen 2003; Shen & Lou 2004a; Lou & Wu 2005; Lou & Bai 2006). Meanwhile, the ubiquitous presence 
of galactic magnetic fields requires a serious and systematic consideration of the magnetized gas disc in a typical disc spiral 
galaxy in order to trace the dynamic evolution in a more realistic manner. In our idealized MHD model, the more massive dark 
matter halo is presumed to be axisymmetric and compatible with the scale-free conditions, while the magnetic field is taken 
to be of an isopedic geometry. Our MHD disc model does include several essential elements of disc galaxies and have already 
shown some interesting and intricate phenomena. In short, self-gravity and thermal pressure of a rotating scale- free disc, the 
axisymmetric gravitational potential of a dark matter halo and the isopedic magnetic field across the disc are synthesized 
together in our current model framework. 

In the model analysis, we consider the self-gravity component in the disc without the usual WKBJ or tight-winding 
approximation and allow the background disc scale-free index /3 and the perturbation scale-free index /3i to be difi'erent in 
general. In other words, the long-range gravitational effect has been taken into account in full. Therefore equation (1501 is fairly 
comprehensive and rigorous linear dispersion relationship, from which much information can be gained, such as the WKBJ 
relationship in Appendix B, global stationary MHD dispersion relations l|56^ and l|58|l . By choosing different parameters, 
equation H5()|l is sufficiently general to include aligned {v = 0) and unaligned 7^ 0) cases, as well as axisymmetric (m = 0) 
and nonaxisymmetric (m 7^ 0) cases. 

The MHD discs are specified in a range for the background scale-free index j3 £ (—1/2, 1/2) [n.b., /3 £ (—1/4, 1/2) for 
warm discs only], the inverse- square magnetosonic Mach number (disc temperature and magnetic pressure together) w > 0, 
the ratio of the gravity force of the dark matter halo to the unperturbed disc gravity force / > and the gravity dilution 
(due to the magnetic tension force) factor e ^ 1. Meanwhile, MHD perturbations are characterized by the perturbation 
scale-free index /3i £ (—1/2, 1/2), azimuthal wavenumber m (integer) and radial wavenumber i/. With all these parameters 
in physical ranges, we proceed to construct global stationary MHD perturbation configurations in a specific background disc 
(characterized by /3, u), / and e parameters) under a specific perturbation (characterized by /3i, m and v). In comparison 
with the WKBJ dispersion relation (Appendix B), the marginal instability of axisymmetric MHD disturbances is examined 
for discs with fiat rotation curves {l3 — 0). 

Our conclusions are summarized below. 

(i) Aligned Cases with v = 0. 

For global stationary aligned MHD perturbations, we take /3i — [3, although the disc system does allow for I3\ ^ {3 
in general. The axisymmetric (m = 0) aligned MHD perturbation represents purely a slight rescaling of the axisymmetric 
background (Lou & Shen 2003; Shen & Lou 2004a; Lou & Zou 2004, 2006; Lou & Wu 2005). For both e > and e < regimes, 
the global m ^ 2 stationary MHD perturbation configurations exist when w varies within a certain limited region (see Figs. 
1, 2 and CI). For neutral m = 1 MHD perturbation modes with e > 0, there are three types of relations between w and / 
when P parameter falls within the three intervals (—1/2, —0.2071), (—0.2071, 0) and (0, 1/2), respectively (see Figs, la, 2a, 
2b and Fig. C2). When e < 0, neutral m = 1 MHD perturbation modes can only exist for 13 £ (—1/2, —0.2071) (see Figs, lb 
and C2). Having specified values of m, (3 and e parameters, the global stationary MHD perturbation configurations can be 
displayed in the (w, /) diagram. One point in the {w, f) diagram then corresponds to one specific global stationary MHD 
perturbation configuration, and all of these points together constitute continuously a hyperbolic curve in the {w, /) diagram. 
A detailed description and figures can be found in Appendix C. 

(ii) Unaltgned Cases with 1/ ^ 0. 

To maintain a constant radial flux of angular momentum associated with MHD perturbations (e.g., Goldreich & Tremaine 
1979; Shen et al. 2005), the perturbation scale-free index [3i should be set to 1/4. Furthermore, the /3i = 1/4 stationary MHD 
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perturbation configurations, in which the angular momentum flux is constant (see Appendix G), restrict parameter w to a 
certain limited region for m ^ 1 neutral MHD perturbations. In particular, the m = 1 neutral MHD perturbation mode 
does not exist for e < 0. For given values of m, /3, v and e parameters, global stationary MHD perturbation configurations 
correspond to a hyperbolic curve shown in the (w, f) diagram (see Fig. 3); this is similar to the aligned MHD perturbation 
configurations. Details of analysis can be found in Appendix D. When /3i 7^ 1/4 (i.e., non-constant radial flux of angular 
momentum), the dispersion relation becomes complex and its imaginary part gives an additional constraint. Thus for given 
values of m, P, Pi, v and e parameters, the continuous point set of global stationary MHD perturbation conflgurations shrink 
from an extended one-dimensional curve to a single point in the (w, f) diagram (see Figs, la, 2a, 2b). Numerical explorations 
show that only m = 1 lopsided MHD perturbation modes exist. 

(iii) Marginal Instabilities for Axisymmetric Disturbances with m = Q, /3i = 1/4 and v ^ 0. 

Given a special value of (3 — 0, we have shown that the increase of a dark matter halo mass (larger / parameter) or of 
magnetic field flux (smaller e > parameter) can squeeze both the Jeans collapse and ring fragmentation instability regions. 
In the regime of e < 0, the MHD disc becomes stable for axisymmetric MHD perturbations of all radial wavelengths, according 
to the WKBJ analysis. 

From the model results summarized above, we can see the existence of an isopedic magnetic field can greatly influence the 
global stationary MHD perturbation configurations in a scale- free disc system, as well as the axisymmetric stability properties. 
When < e < 1, an isopedic magnetic field gives rise to an effective modification factor on the gravitational constant G, from 
G to eG (Shu et al. 2000), and an effective enhancement factor on the sound speed Vs (thermal pressure), from Va to \/Qvs. In 
the absence of the dark matter halo, e has to be positive in order to maintain a disc equilibrium. In other words, the magnetic 
field tension force cannot dominate the self-gravity component in the disc plane. However, a massive dark matter halo exerts 
an extra gravity force to hold the magnetized disc together, allowing the magnetic tension force to overtake the self-gravity 
component in the disc plane, corresponding to a situation of e < 0. In the parameter regime e < 0, we can no longer regard 
the isopedic magnetic field as a modification but must face an entirely new physical situation involving several interesting 
features in a strongly magnetized disc system. An interesting example is to derive the MHD WKBJ relation in equation l|52^ 
for the instability of axisymmetric disturbances (m = 0). When e < 0, we always have uP' > 0, indicating the disc system 
being stable for axisymmetric MHD disturbances of all wavelengths. 

Starting from the case of e < 0, we may extend the same idea to more generalized situations. As a system has a strong 
magnetic field attached to plasma or gas, the magnetic field would have torn the system apart in the absence of an external 
massive dark matter halo potential. Two possible physical examples come into mind: the hot magnetized gas trapped in the 
centre of a galaxy cluster (e.g., Makino 1997; Dolag et al. 2001; Hu & Lou 2004) and the accretion disc around a SMBH 
(e.g., Kudoh et al. 1996; Kaburaki 2000) where strong magnetic field may exist. The former can be perceived in reference 
to a quasi-spherical system in which a hot gas is mainly confined by the dark matter potential and the hot gas itself is 
strongly magnetized. Our current model offers another possibility in that a more or less fiat gas disc system in rotation can be 
primarily confined by a dark matter potential and the gas disc itself is isopedically magnetized with a considerable strength. 
By observing the existence of unusually strong magnetic field in a localized region, we may predict some 'unseen' dark matter 
components of the system. For instance, if a source region has strong cyclotron or synchrotron emissions with much reduced 
optical emissions, the scenario outlined above may provide a possible explanation. 

Fujita & Kato (2005) proposed that Weibel instability may be responsible for generating strong magnetic fields in galaxies 
and clusters of galaxies at redshift as high as 2 ~ 10; in such a scenario, the magnetic field can play a crucial role in forming 
galaxies and galaxy clusters. Our model results, especially those for a strong magnetic field (i.e., e < 0), offer several interesting 
clues for the dynamics of galaxy or galaxy cluster formation. With extensive observations for magnetic fields in clusters of 
galaxies (Taylor et al. 2006), a more sensible MHD model may be constructed and tested. 

The final point of interest involves the radial flux of angular momentum. If the radial flux of angular momentum is 
constant, only the central part will collapse and the outer portion of a disc is stable. If the radial flux of angular momentum 
increases with increasing R, the disc mass at all radii will lose their angular momentum by interacting with stationary MHD 
density waves, and hence gradually spiral inward and drift towards the centre. In this situation, not only the central part, 
but also the entire disc drifts inward. For a single rotating disc embedded in an axisymmetric dark matter halo, the angular 
momentum of the disc system is conserved. In the presence of an external nonaxisymmetric potential, such as a companion 
or satellite galaxy, this conservation no longer holds. Imphcations of this issue will be investigated further. 
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APPENDIX A: SEVERAL PROPERTIES OF FUNCTION YMiP') 
The function Ym(/9') can be explicitly expressed as 

' r(2i + i/2)r(z2 + 1/2) 

in terms of Gamma function r(- • ■ ), where 

zi = m/2 - + 1/2 , 22 = m/2 + P' , 13' = - iv/2 . 

Since T{z*) — T{z)* with the asterisk * denoting the complex conjugate operation, it then follows that 

7vr{zi)r{z^) _ 7rr(2i)*r(22)* 



7 = Yr^iP'T 



r{zl + l/2)r(22* + 1/2) T{zi + 1/2)T(22 + 1/2) 
Using the recursion relation r(2)r(l — z) — 7r/siri(7rz), we obtain 
Y_^{l3') = ^r(l/2-22)r(l/2-2:i) ^ YraiP') ^^['^(^ ^ ^2^] sin[7r(l - zi) 



r(i-z2)r(i-zi) 



.[7r(l/2-Z2)]sin[7r(l/2-2i)] 



cos TT I 21 — 22 ) 



COS 7r(2l + 22) 



That is, Ym{P') function is symmetric or even with respect to the subscript m. 

According to the Stirling formula, the asymptotic series of Gamma function r(2) is given by 

r(2) = {2n)'^^e-'z'-^^'' (1 + — + + ■■■] . 

^ ^ ^ ^ V 122 2882^ J 

In terms of asymptotic series expansion 1A4II . r(2i), r(22), r(2i + 1/2) and r(22 -|- 1/2) are given by 



r(^i) 


= (27r)'''^e 


r(22) 


= (27r)'/''e 


r(2i + 1/2) 


= (27r)'/''e 


r{z2 + 1/2) 


= (27r)'''^e 



21-1/2 

^1 



1 + ^+0(0 

i + ^ + o(.-) 

To the leading order, we then have Ym{P') in the form of 
y„(/3') = ETTzt'-'^^z'^'-'^^izx + l/2)-^i (22 + 1/2)'^^ [1 + 0{u-^)] . 
A formula frequently used in the following analysis is given below 
z' = exp [2'(ln I2I + i arg 2)] , 

where 2 and 2' are two arbitrary complex numbers and arg 2 is the principal value of the argument 2. 
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In addition, we have 



ln|z, + 1/2| 
In I Zi 

and 

argzi 
arg 22 
arg(2i + 1/2) 
arg {Z2 + 1/2) 



in I V'[3ff(^> + l/2)]2 + [9(2, + 1/2)]2| =\n(\v\/2) + 2[^[z, + l/2)f/v'^ + 0{v-^) , (fori = 1,2) 
in {ypPOFTpPOF} = HH/'2) + 2mz,)f/i^^ + 0(1^-") , (for i = 1, 2) 



sgn(//)7r/2 - 2U(zi)/i' + 0(u-^) , 
= -sgn(!^)7r/2 + 2K(Z2)/!^ + 0{u-^) , 

sgn{u)-K/2-2^{zi + l/2)/iy + 0{v~^) 



-sgn(i.)7r/2 + 2K(22 + 1/2)/// + 0(i/"^) , 
where sgn stands for the signum function. Hence, we arrive at the following two relations 

Li = (21 - l/2)lnl2il + (22 - l/2)ln|22l - 2iarglnl2i + 1/2] - 22hil22 + 1/2] = - ln(li/l/2) - i(l - 4/3i)/2// + 0(1/"^) , 
and 

L2 = (21 - 1/2) arg 21 + (22 - 1/2) arg 22 - 21 arg (21 + 1/2) - 22 arg (22 + 1/2) = i + (1 - 4/3i)/i/ + 0(z/"^) . 

From formula 1A6II . it is easy to see 
exp(Li + iL2) = zl'-^^'^z'^'-^^^z^ + l/2y (22 + 1/2)""^ . 

Finally, we derive 



2tt 



i(l-4/3i) 



2^ 



(A7) 



APPENDIX B: THE WKBJ OR TIGHT- WINDING APPROXIMATION 

According to equation we have 

(^ + ^fe„)^_2(l-/3)6g ^ ,^+^^+2/3.-4A^+iJ-- - 2(l + /3)^(^ + mM 



wefeg -ecoGKn(/3') 



4/3i - 1 + 



2/3mbo 
Qj + mfeo 



which can be further cast into the form of 
i?2[(t^ + m!r2o)^ 



^ 

ew2g-ei?EoGY„(/3') 



= + in + 2/3i - 4/3i + i;/ 



(ii + m6o)2-2(l-/3)fo2 

4(1 + /3) [T7ibo - /3i (1:; + mbp)] ci) 
(ii + m6o)2-2(l-/3)fo2 

2(1 + /3)a;(aj + mSlo) 



4/3i - 1 + 



(oj + mfio)^ 



+ 



2/3mf2o 4(1 + /3) [mfio - /3i (a; + mfio)] 



(a; + mfio)^ — Kg 



(Bl) 



^.2 



0. These correspond to the corotation 



First, it is clear that equation (IBH is singular when uj^mVL = and (ctJ + mf^o) 
resonance and the inner and outer Lindblad resonances (Goldreich & Tremaine 1979). In our analysis, such singular points 
do not appear due to the stationary condition lj = 0. 

In the regime oi v max(m, c^;, 1), equation IBlll can be significantly simplified to 



(a; + raVL^Y - Kq 
ew2g-ei?EoG'y„(/3') 



R 



[1 + 0(^.-^)1 



(B2) 



From Appendix A, we know ym(/3') ~ 2-k l\v\ in the limit of ^ max(m, 1). By defining the radial wavenumber k as v jR, 
we can rearrange equation 1B21 into the following form 



(oj + mS7o)^ = kq + fc^Ofso ^ 27reG|fc|Eo 



(B3) 



which shares the classic WKBJ form (Lou & Fan 1998a). According to equation l|B2|l . dispersion relation ljB3^ has an error 
of the order of 0{y~^^. From equation IBIL we note that errors can be reduced to 0{y~^^ when /3i — 1/4 and lj ^ 0; for 
these special cases, a more accurate dispersion relation can be obtained (Shu et al. 2000). 
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f-^{-^-C2] , (Cl) 



APPENDIX C: ALIGNED CASES 

For aligned cases, we begin with the expression 

W — Ca 

where the three coefficients Ca, Ci, and C2 are exphcitly defined by 
^ m^-2(l-/?) 

- m2+ 4/3(1-/3) ' ^^^> 
_ [l + (l + 2/3)a]y„(/3) 

= — mm — ' ^^'^ 

_ ^ (l + 2/3)y^(/3) 
= ' 2/3Fo(/3) ■ 

Clearly, we can write equation IClll in the form of (/ + eC2){u> — Ca) = eCi. When plotted on a two-dimensional diagram of 
(w, /), this relation takes a hyperbolic shape. 

The three constraints w ^ 0, f ^ and / + e can also be effectively written as ^ and / ^ max{0, — e}. 

In short, equation ICH depicts a hyperbolic curve in the {w, f) diagram for specified m, (3 and e parameters. The two 
constraints w ^ and / ^ max{0, — e} enclose a region of physical interest. In other words, only the enclosed parts of the 
hyperbolic curve are physically meaningful. 

For different combinations of m, (3 and e, we will discuss whether the physically meaningful solutions of the hyperbolic 
curves exist and the possible ranges of the curves, which are expressed as ranges of w. From the hyperbolic curves, the ranges 
of / can be readily calculated for the given ranges of w. Note that m is a positive integer, /3 G (—1/2, 1/2) and e ^ 1. The 
special case of e = is investigated in the main text. 



Cl Cases of m ^ 2 

For cases of m ^ 2, our numerical calculations show that 1/2 ^ Ca < 1, Ci > and 1/3 < C2 < 1. The case Ca ~ 1/2 is 
achieved for m = 2 and /? = 0. 

When < e ^5 1, the constraints are / ^ and w ^ [see panel (a) of Fig. Cl]. We can further infer 



/ > , ™ > ^ , , >C2 W e (Ca , Ca + ^ 

When e < 0, the constraints are / + e ^ and w ^ [see panel (b) of Fig. Cl]. We can further infer 



(C5) 



/ + O0, w^O max <j Ca + ^^^^ , 0} ^w<Ca w e [0, Ca) . (C6) 



C2 Cases of m — 1 

When m = 1, both coefficients Ca and Ci may become sing ular at /3 = (1 - a/2)/2 = -0.2071. It follows that C„ > for /3 G 
(-1/2, -0.2071) and Ca < for /3 G (-0.2071, 1/2), respectively. We also have Ci > for /3 G (-1/2, -0.2071) U (0, -1/2) 
and Cl < for /3 G (-0.2071, 0), respectively. Finally, C2 < for /3 G (-1/2, 0) and C2 > for /3 G (0, 1/2), respectively. 
We now consider several subcases separately. 



C2.1 The subcase of 13 € (-1/2, -0.2071) with Ca > , Ci > , C2 < . 
When < e s; 1, 

/ > , w^O ^ wG{Ca, +00) , (C7) 

and the range of w can be identified in panel (a) of Fig. C2. 
When e < 0, 

f + e^O, w^fO max |ca + ^^^ZJ > o| sC ^ Ca ^ w G [ 0, Ca ) , (C8) 

and the range of w can be identified in panel (b) of Fig. C2. 
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Figure CI. Coiiibiiiatioiis of w (a measure for the disc temperature and magnetic pressure together) and / (a measure for the background 
dark matter halo potential strength) for m ^ 2 global stationary MHD perturbation configurations when /3 € (—1/2, 1/2). Panel (a) 
shows the case of e > (i.e., the disc self-gravity overtaJses the magnetic tension force) and panel (b) is the case of e < (i.e., the 
magnetic tension force overtakes the disc self-gravity). The physically meaningful parts are plotted in solid curve and the unphysical 
parts are plotted in dashed lines. In panel (b), the dash-dotted line indicates the constraint on / when e < 0, due to / -|- e ^ 0. Several 
asymptotes are given in dotted lines, / = —eC2 and w = Ca- 



C2.2 The subcase of P € (-0.2071, 0) with Ca < , Ci < , C2 < . 
When < e ^ 1, we infer 

/^O, w^O ^ w ^ max|c„-|- ^ , o| ^ w e C„ -h ^ , -|-oo^ , (C9) 

and the range of w can be identified from panel (c) of Fig. C2. 
When e < 0, wo infer 

2 

f + e ^ , tZ)>0 => Ca ^ w ^ Ca + -p^ r = ^ , w > => Such a w cannot exist , (CIO) 

G2 — 1 1-1-2/5 

and no curve is within the physically allowed region in panel (d) of Fig. C2. 
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C2.3 The subcase of 13 £ (0, 1/2) with Ca < , Ci > , C2 > 
When < e ^ 1, we infer 



/ > , ^Z; > 



max{Ca , 0} ^ < Ca + 



Ci 
C2 



w G 



0, c„ + 



C2 



and the range of w can be seen from panel (e) of Figure C2. 
When e < 0, we infer 

/ + e^0, w ^ =J> i(;<Ca, w ^ Such a w cannot exist , 

and no curve is within the physically permitted region in panel (f) of Figure C2. 



(CU) 



(C12) 



APPENDIX D: SEVERAL PROPERTIES OF /3i = 1/4 DISCS 

For a disc system with /3i — 1/4, the global stationary dispersion relationship for MHD perturbations can be transformed 
into 



where 
Do = 

Di = 

Do = 



[m" - 2(1 - m 
(m2+2/3 + l/4 + i^2) 



w ~ Do 



-A 



(/ + eD2)(™-Do) = eA , 



2/3^0 



r(m/2 + l/4 + ii//2) 



r(m/2 + 3/4 + i^/2) 



[2m2(l + /3)+ 4/32 -7/4 + 1/2] 



(m2 + 2/3 + 1/4 + 1/2 



7r(l + 2/3) 
2/3Fo(/3) 



r(m/2 + l/4 + ii//2) 



r(m/2 + 3/4 + ii//2) 



(Dl) 

(D2) 
(D3) 
(D4) 



Numerical computations show the following results: —8 < Do < for m = 1 and < Do < 1 for m ^ 2, respectively; in all 
circumstances, D\ ^ and Di = occurs when m = 1, /3 = —1/4 and i/ = 0; in all circumstances, ^ D2 < 1 and D2 = 
occurs when m = 1, /3 = -1/4 and i/ = 0. Note that /3 £ (-1/2, 1/2). 

Given m, v, (5 and e parameters, equation HDH represents a hyperbolic curve in the {w, f) diagram. Two physical 
constraints are w ^ and / J? max{0, — e}, enclosing the physically meaningful parts of the hyperbolic curve expressed as an 
area of w. 

The case of e = is special, in which / becomes arbitrary and w — Dq. As constrained by w 0, neutral m = 1 MHD 
perturbation mode cannot exist in e = discs. 



Dl Cases with m ^ 2 

In these cases, we have < Do < 1, Di > and < D2 < 1 for m ^ 2. When < e ^ 1, the two physical constraints are 
/ ^ and w ^ 0. We therefore infer 



f>0 



Dl 



^ D2 



w - Do 

When e < 0, the two physical constraints are / + e ^ and w ^ 0. We then infer 

Dl 



w) G Do , Do + 



Dl 
D2 



/ + e>0 



w>0 



max < Do 



0} i^w < Do 



we[0, Do) 



D2-1 

The schemata of the above two cases (e > and e < 0) are the same as panels (a) and (b) of Fig. CI, respectively. 



(D5) 



(D6) 



D2 Cases with m — 1 



In these cases, we have Do < 0, Di > and < D2 < 1 for m = 1; the Di = and D2 = cases will be discussed separately 
below. When < e ^ 1, we infer 



/ ^ , w^O 



Dl 



> D2 



w € 



, Do + 



Dl 
D2 



w — Do 
When e < 0, we infer 

f + e ^ , w ^ w < Do , w ^ Such a w cannot exist . 

The schemata of the above two cases (e > and e < 0) are the same as in panels (e) and (f) of Fig. C2, respectively. 



(D7) 
(D8) 
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APPENDIX E: LOGARITHMIC SPIRAL DISCS WITH u^O AND /3i / 1/4 

When u ^ and /3i 7^ 1/4, a point in the (it), /) diagram is given by the two values w and / explicitly expressed as 

where 5R(- ■ • ) and • ■ ) represent the real and imaginary parts of the argument, respectively, 

. . ^ [m^ - 2(1 - /3)](m^ + 2/3 + 2ft - 4/3? + /.^) 
^ ~ (m2 + 2/3 + 2/3i- 4/32 + 1.2)2 + j,2(i_4^^)2 ' 

c>^n = Kl"4/3i)[m^- 2(1-/3)] 

~ (m2 + 2/3 + 2/3i- 4/32 + 1/2)2 + j,2(i_4^^)2 ' 

In the limit of ^ +00, we have the following asymptotic expressions to the leading orders according to equations iA7l . iE4l 
and 1E5I I. namely 

WmiP')] = ^[1 + 0{U'^)], (E6) 



c,[YM - ^ii^^lM[i + o(^-)], (E7) 



^3 

K(C) = !I^l^l(l^[l + 0(^-^)] , (E8) 
^iC) = (^^^^-)t-r^(^'^)][l + 0(.--)]. (E9) 



To the leading order of large i', we then have w and / expressed as 



(ElO) 



- - 2/3yo(/3)[:2j2(l-/3)] [^ + ^(-")i- 

According to equations l|E10|l and l|Ell|l with 1/ — > +00, w < when m ^ 2 and w > when m = 1, respectively. For m = 1, 
/ > when e > and / < when e < 0, respectively. 

The above analytical analysis based on the 1/ — > +cx approximation gives the variation trends of w and /. But, this 
approximation becomes invalid when ~ 1 or smaller. We resort to numerical experiments to probe solution properties. 
According to extensive numerical tests, we find that 5R[ym(/3')] > and !/(l — 4/3i)Sj[y,„(/3')] > when m ^ 1. 

El Cases with m ^ 2 

When m ^ 2 , we have 5R(C) > . Meanwhile, w can be written as 

im^ + 2p + 2/3i - 4/32 + .2) c,[y„(^,)] ] ■ (E12) 
Extensive numerical results show that 

Kl-4/3i) g?[y^(/3')] Kl-4/3i) , 

(m2 + 2/3 + 2/3i - 4/32 + i/2) SJ[y,„(/3')] (m2 + 1 + 2/3i - 4/32 _^ ^,2) 9[y;„(/3')] ^ ' 

Therefore when m ^ 2, we have < by extensive numerical explorations. 

E2 Cases with m = 1 

In these cases, we have 

- / 1 + 2/3 + 2/3, - 4/3? + - .(1 - 4/3,) 



[(1 + 2/3 + 2/3i- 4/32 + 1.2)2 4. j,2(i _ 4^^)2] ' ' ^i' v c^fy^ (/j')] 

Extensive numerical explorations show the following inequalities 

1 + 2/3 + 2/3, - 4/3? + - .(1 - 4/3i)|jm;j| < 2 + 2/3, - 4/3? + - .{I m ^^^^l^^l < 
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I'D 1.123 1.761 



Figure Fl. The axisymmetric (m = 0) marginal instability curves (upper-right and lower-left solid curves) of i^/w) (the reciprocal 
of the disc temperature and magnetic pressure together) versus u (dimensionless radial wavenumber for perturbations) for specified 
values of f/e = 0.5, the scale-free index /3i = 1/4 for perturbations and the background scale-free index /3 = 0. Here, e is taken to be 
positive. Five specific points, uq, 1.123, 1.760, Us and ui,, can be found along the p-axis. The dashed line of Us gives the asymptote and 
the dash-dotted line marks the inflection point i/f,. 



and therefore we have w > . 

The / parameter can be written as 



[1 + (1 + 2/3)iD]K(B) 



- 1 



[w ~ 

We can also derive the following relationship 

3fJ(B)Si(C) _ !/(l-4/?i)(l-2/?) 



(E13) 



^ ' 9(B) (1 + 2/3 + 2/31-4/3?+ 1.2)2 -f 1.2(1 _4^^)2cj[y^(^,)] 

For e < 0, we require / + e ^ on the ground of physics. However, 
[1 + (1 + 2/3)w]iR(B) ) 



> 



< 



K(c)] 7 

leading to an obvious contradiction. Therefore for e < 0, stationary m = 1 MHD perturbation modes do not exist. 
When e > 0, the physical constraint becomes / S5 0, giving the following inequalities 



[1 + (1 + 2/3)u;]5R(B) 



^ 1 



(1 + 2/3)5R(C) + 1 



S5(C) 



(l + 2/3)5R(B) - 1] > 



Extensive numerical examples show that the above inequalities always hold true. We therefore infer numerically that when 
e > and / > 0, stationary m = 1 MHD perturbation neutral modes do exist. 



APPENDIX F: AXISYMMETRIC CASES 

When m = 0, /3i = 1/4 and /3 = 0, we have 



w ~ AAo(i.)-8£:/(l + 4i.2) ' ^ 

where the Kalnajs function Noi^v) is given by definition H15|l and £ is defined as 

£ = l + f/e. (F2) 

Since we require inequality / + e > to hold, it follows that f = l + //e<Oife<0. When £ <Q, w remains always negative 
and thus unphysical. For this reason, we only focus on the situation of e > and thus £ ^ 1. 

When £—Mo(v) — 0, one has = by equation ijFl^ . Numerical explorations show that No{v) decreases with increasing 
V. When £ > No(Q) = 4.377, a positive £ — Noif) can be always guaranteed. We adopt the notation uq for the specific v 
making £ — NoIuq) — 0. Hence Afo{i'o) — £ ^ I corresponds to i/q ^ 1.123 because ^0(1. 123) = 1. When £ > 4.377, there is 
no such uq root. 
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The singular point Vs is to make Ao (i^s ) — 8.5 / ( 1 + ii^g ) = in equation HFH . Extensive numerical explorations show that 
increases with increasing v. Therefore, the singular point fs should be larger than 1.760 which is the root of 
equation Afoii^) - 8/(1 + Au^) = 0. 

From equation IFlfl . we readily obtain 

dD^ _ (l + 4i/^)AAo(t^)(4i/^-7) 

Thus as £ = 1 + //e increases, grows for u > \/7 /2 — 1.323 and decreases for < i/ < \fl /2 = 1.323, respectively. 
Meanwhile, we can also readily derive 

dj^ ^ Uv£ L M - AA(;M(l + 4t.^)(4t.^-7) _ ■ 

We define a special value ui, which is the solution u of the following equation 

AAoM-^^iMilMM^^.^O. (F3) 

Numerical explorations clearly indicate that Ub ^ 3.058 which is the root of equation <F3t with £ — 1. With the increase of 
u, decreases when v < Vb and grows when u > Ub, respectively. Moreover, numerical results show Ub > Us- 

Finally, we have an order of inequalities vo ^ 1.123 < 1.323 < 1.760 ^ Vs < i^b (see Fig. Fl). When i/ G [0, vo] [J{j^s, +oo), 
we have ^ 0. More specifically, is a decreasing function of f when u £ [0, i^o] UC'^s, Vb) and an increasing function of 
v when u G (ub,+oo). In addition, is a decreasing function of £ when u G [0, uo] and an increasing function of £ when 
I' G (i^s, +oo). When £ > 4.377, both uq and the range [0, i^o] disappear. Since £ = 1 + //e, the variation trend of as / 
and/or e change can be readily determined. 



APPENDIX G: ANGULAR MOMENTUM FLUX TRANSPORT 



In this appendix, we provide several proofs and discussions for the angular momentum flux transport associated with global 
coplanar MHD perturbations in a scale-free disc system with an isopedic magnetic field geometry. In general, such angular 
momentum flux transport contains three separate contributions (Lynden-Bell & Kalnajs 1972; Fan & Lou 1999; Shen, Liu & 
Lou 2005), namely, the flow advection transport deflned by 

f2TI 



the gravity torque flux transport A*' defined by 

.G_ 1 



A^ 



AttG 



d9 



dzR 



do 



9^(01) 



dR 



and the magnetic torque fiux transport A defined by 



de 



d2K(6fl)K(be) 



(Gl) 



(G2) 



(G3) 



The three-dimensional gravitational potential perturbation (j)\{R,9,z) is associated with a Fourier harmonic component of a 
coplanar logarithmic spiral perturbation in the surface mass density Ei = c\R~^^ ~^ exp(im&) = ci 7?"^''"^"^"' exp(im&) = 
C\R~'^^~^ exp[i(m6'-|-!/ln R)] where ci is a small constant amplitude coefficient. By the convention of our notations, <f>i{R, 6*, 0) = 
$1 with "I>i being introduced in equation 14111 . Here, 6n and be represent the two components of the magnetic field perturbation 
in the isopedic magnetic field. From equation (2.14) of Shu & Li (1997), we have 

1 901 ,1 d(j>i 



By definition 27rGA = \G^^^ and l/A^ 



A'= 



dz-R{bR)-R{be 



2-kGA dR ' 
1 — e, we then have 
1 



47rGA2 



de 



2-kGKR de ■ 



dzR 



de 



aaf?((^i 

OR 



(e-l)A^ 



(G4) 



(G5) 



Corresponding to Ei = ciR '^^ ^ exp[i(m& + u\nR)], the gravitational potential perturbation (t>\ can be written as (e.g. 
Binney & Tremaine 1987) 



i(i?,6',z) = -27rGciexp(im6i) / dkeyi-p{-k\z\)J^{kR) dR R 



'"JmikR') , 







where Jm{x) is the cylindrical Bessel function of order m with an argument x. We then have 

6i)=-27rGci/ dkexp{-k\z\)J^{kR) dR' R'''^^^ cos{me + v\nR')J^{kR') 
Jo Jo 



(G6) 



(G7) 
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where ci is taken to be real. 
It follows immediately that 
a3fJ($i) 



de 



2-KGmCl I dkaexp{ — ka\z\)Jm{kaR) I dRaRa^'^^ Sm{m9 + V In Ra)JmikaRa) 

Jo Jo 



(G8) 



and 

^^^A ^ dkte^p{~kb\z\)\^^kbJ„,+i{kbR) + ^JmikbR)\ J^^ dRtR^""^' cos{me + uln Rt)Jm{kiRb) . (G9) 

Using the following relations 

2 



/OO P OO 

dzexp[— (fca + kb)\z\] = 2 / dzcxp[— (fca + kb)z] = — — 
-OO Jo (ka 



and 



_ 1 

^ 2 

we obtain the gravity torque flux transport 



/ d&cos(mS + In /?(,) sin(m& + i^ln 
Jo 

d6l{sin[z/ln(_Ra/-R(,)] + sinprn^ + v\n{RaRb)\} = vrsin 



V In 



(GIO) 



(Gil) 



2m-K^GciR I dkadkbdRadRbR-^^^ R'^^'^KJrr.(kaRa)Jm{kbRb) 



V In 



Jm{kaR ) \m 

{ka + kb 



Jm(fcb-R) - kbJm + l(kbR) 



By taking the following integral transformations 



Xa = kaR , Xb = kbR , Va = ^ , TJb = ^ 



ka = ^ , kb = ^ , Ra = VaR , Rb = VbR 



we arrive at 

A'^ = -~2rmT^GclR^-^'^^M{f3i,m,u) 
where the function m, is defined by 



Mi(3i,m,u) = / dXadxbdyadybivaVb)'^'^^ Jniixaya)Jm{xbyb) sin 
Jo 



V In 



(G12) 

(G13) 
(G14) 

[mJm{xb) ~ XbJm + liXb)] ■ (G15) 



Xa ~t~ Xb 

We correct here a few errors in Appendix C of Shen et al. (2005): Jm{k'r) is redundant in their equations (CI), (C2) and 
(C4); Jm{x') is redundant in equation (C6) and in the expression of Ep below equation (C8), and the heavy solid curve of H*"" 
in their Figure CI is thus invalid; an extra numerical factor 2 should be multiplied in their equations (C4), (C6) and (C8), 
respectively. 

According to equations 16011 and 1611 1. we have 

= 2mcii?~^''^'''''/C_Riexp[i(m^ + vlnR)] , vei = —ciR'^'^^^'^ICe exp[i(m6' + vlnR)] , 

where 



(G16) 



ICr = 



[wObf, ~ ecoGYmiP')] (1-/3') 



ICo 



[wQbl - ecoGYmjP')] [m^ - 2/3' (1 - /3)] 
cofoo [m2 - 2(1 - /3)] 



cobo [m2 - 2(1 - /?)] 
We therefore come to 

?lt{vRi) = ~2mciR~^'^^+'^ 1^{ICr) sin{m6 + iy\nR)+ Q{JCr) cos{me + ulnR)] , 

and 

3?(«ei) = -cii?-^'^i+'^ [?St{K:e)cos{me + ulnR) -Q{K:e)sm{me + ulnR)] . 

Since 

/ de [K(/Cjj) sin(m6' + i^lnR) + Q{ICr) cos{me + i^ln R)] [K(/Ce) cos(m6' + In i?) - 9(/Ce) sin(m.6' + In R)] 
Jo 

= o / {-SR(/CH)9(/Ce) [1 - cos(2m6' + 2iy In R)] + 3ff(/C9)9(/CH) [1 + cos(2m6' + 2iy In R)]} 
2 Jo 

= TT [K(/Ce)9(/CH) - K(/Cfl)9(/Ce)] , 

we derive in a straightforward manner that 

A^ = i?^Eo / de^{vRi)U(vei) = 2mnclcoR^~*'^' [U{K.e)Q{K:R) - U{K:R)Q{K:e)] . 



(G17) 
(G18) 



(G19) 
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Finally, we write the total angular momentum flux Atotai as 

Atotai = A* + A^' + A^ = A^*" + eA'^ = 2mR^'*^^cli^Tvco [K(/C9)9(/Ch) - ?fiiICR)QiK.g)] - en'-'GMiPi, m, u)^ . (G20) 

As expected when m = 0, we have Atotai = 0. When (3i = 1/4, Atotai is independent of R. We note that all these foregoing 
analyses axe performed under the stationary assumption of a; = 0. 

This paper has been typeset from a T^JC/ MJjX file prepared by the author. 



